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ABSTRACT 



Procedures to Predict Vertical Differential Soil Movement 
For Expansive Soils. (December 1997) 

Donald David Naiser, Jr., B.S., Texas A&M University 
Chair of Advisory Committee: Dr. Robert L. Lytton 

Damage to lightly loaded structures, paving and service piping in areas of 
expansive clay soils has occurred throughout the world. The cause of this damage has been 
the inability to accurately model expansive soil movement so that foundations are 
adequately designed to withstand the movement. The amount and rate of differential soil 
movement for expansive soils is due to a combination of soil characteristics, namely: 
suction compression index, unsaturated permeability, and diffusivity. Currently, 
geotechnical engineers run tests to measure the soil properties required to estimate 
differential soil movements. However, there seems to be apprehension toward attempting 
these soil movement calculations due to the perceived complexity of the calculations or a 
simple lack of understanding of the theory. The procedures delineating the step by step 
process used to calculate suction profiles and volume strains of expansive soils is 
presented. These procedures include the methodology to predict soil heave and shrink 
underneath shallow foundations which generate maximum center lift and maximum edge 
lift slab distortion modes. The main contributions of this research are: equations and 
procedures to calculate the equilibrium suction profile and depth to constant suction for a 
particular soil profile and location, equations to calculate the horizontal velocity flow of 
water in unsaturated soils, the methodology to predict differential soil movement shortly 
after a slab has been constructed and before the soil under the slab has reached an 
equilibrium moisture content, and the procedures to apply differential soil movement 
theory to soil profiles with multiple layers and moisture effect cases to be used for shallow 
foundation design. 
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INTRODUCTION 



When a slab barrier is placed on an expansive soil, the immediate effect is that 
evaporation and precipitation cannot effect the soil beneath the slab. During subsequent 
weather events occurring to the surrounding soil, moisture differentials develop between 
the uncovered soil and the soil beneath the slab barrier. These differentials are largest 
when the climate is described by a distinct wet period followed by a hot dry period. Unless 
a foundation system is installed to prevent or reduce moisture movement, moisture 
redistribution occurs within the soil mass to reduce these differences in moisture content 
(Holland and Lawrance 1 980). Damage caused by this type of moisture differential is 
considerable and warrants investigation in developing a more scientific and accurate 
methodology to address these problems. 

The ability to predict differential soil movement is critical for engineers to 
formulate a picture of what the soil/structure interface will be. This information is 
necessary to adequately design foundations or other ground supported structures. In slab- 
on-ground design, the principal interest is in making an accurate estimate of the range of 
vertical movement that must be sustained by the foundation. The procedures contained 
herein are methods for defining envelopes of maximum heave and shrinkage so that 
adequate slab-on-ground designs are made. For highway and airport pavements, canals, 
and pipelines, the wave spectrum of differential movements versus wave lengths are the 
desirable design characteristics. Structural floors suspended above expansive clays must 
be provided with a gap that exceeds the total expected heave. Drilled piers (or shafts) must 
be designed to resist simultaneously a vertical movement profile and a horizontal pressure 
profile, both of which change with wetting and drying conditions. Retaining structures, 
basement walls, rip rap, and canal linings must be designed to withstand lateral 
movements. Finally, all foundations must be designed against the time-dependent vertical 
and horizontal curvature that is generated by down hill creep (Lytton 1994). 



This thesis follows the style of the Journal of Geotechnical Engineering. 
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Each type of soil movement listed above and their effects on different types of 
ground structures is of sufficient importance and complexity to warrant research of its own. 
Addressed herein are the effects of vertical differential soil movement on commercial and 
residential types of slabs-on-ground. Vertical differential soil movement in these cases is 
defined as the total amount of vertical movement between two distinct locations separated 
by a characteristic distance. The characteristic distance is described as the edge moisture 
variation distance plus the distance a moisture effect source is away from the slab edge. 

The model used to predict differential movement of expansive soils is based on constitutive 
soil relationships derived primarily from works developed by C. W. Thomth waite (1948), 
W. R. Gardner (1958), R. L. Lytton (1973, 1977, 1994, 1995, 1997), J. R. Pearring (1968), 
J. H. Holt (1969), E. C. Mojekwu (1979), J. L. Nieber (1981), R. G. McKeen (1981, 1992), 
P. W. Mitchell (1980, 1984), D. A. Gay (1994). 

It has been well documented that the magnitude at which soil shrinks and swells is 
not exclusively based on the soils characteristics. The magnitude of differential soil 
movement in expansive soils is dominated by the following variables: soil characteristics, 
pre-construction and post construction vegetation, slope of perimeter surface grade, shape 
and age of slab and depth of soil moisture active zone (Castleberry 1974, Mathewson et al. 
1975, Mathewson et al. 1980). 

The methodology used to predict vertical differential soil movement involves using 
soil suction profiles under the foundation and at the edge of the foundation. These suction 
profiles are developed using the following variables: depth, time, local surface annual 
weather and vegetation conditions, and soil characteristics such as the suction compression 
index, unsaturated permeability, and soil diffusivity. Soil suction profiles combined with 
overburden stresses are used in volume strain equations developed by Lytton (1994) to 
predict the volumetric strain for an incremental layer of soil. The incremental layer 
movements are summed up and total vertical differential soil movement is predicted for a 
particular column of soil. These procedures are detailed in the chapter on “The Procedures, 
in a Software Design Format”. These procedures address foundation perimeter moisture 
effect cases that are either very common with many foundation problems or have a 
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significant effect on vertical differential soil movement. Specifically, the foundation 
perimeter moisture effect cases that are addressed include: (1) bare soil at the surface, (2) 
grass at the surface, (3) flower bed at the surface with a defined moisture effect zone depth, 
(4) trees at the surface with a defined root zone depth, and (5) moisture effect zones that 
have been measured in the form of suction profiles. Included are methods for predicting the 
effects that vertical or horizontal barriers have on vertical differential soil movement. 
Additionally, methods to predict differential soil movement shortly after construction and 
before the soil has reached an equilibrium moisture content under the foundation are the 
highlights of this research and are presented. 

In “The Problem” chapter, the theory and methodology used to predict differential 
soil movement is defined and described in order by which the theory was developed. 

These bits and pieces of theory are applied in step by step processes described in the 
procedures sections of this document to address all of the different moisture effect cases, 
with and without barriers imposed. 
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THE PROBLEM 



The Magnitude and Description of the Problem 

Since early 1950, damage to lightly loaded structures, paving and service piping in 
areas of expansive clay soils has occurred throughout the world. Estimates show that one 
in five Americans are affected by this type of damage while only one in ten are affected by 
major floods. Structural damage caused by expansive clay soils in the United States 
annually exceeds that caused by earthquakes, hurricanes, and floods combined (Jones and 
Holtz 1973.) Damage to residential and commercial structures caused by soil movement is 
estimated at 660 million dollars, almost one-third of the total monetary damage received by 
all structures as a result of soil movement. On a state-by-state basis, California and Texas 
account for 35% of all expansive soil losses in the United States. This becomes obvious 
from a map of the location of much of the expansive soils in the United States, see Fig. 1 . 
Much of this soil is located in highly populated regions of both California and Texas. 




■ High 
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Fig. 1. United States Map of Expansive Soils After Wiggins (1976) 
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When it is realized that most residences and light commercial buildings in these two states, 
which together had 430,000 new home starts in 1978, are constructed with slab-on-ground 
foundations, the need to develop a theoretical design approach based on physical and 
natural law for slab-on-ground design in expansive soil areas becomes evident. 

The problem with slabs-on-ground built on expansive soils becomes manifested as 
excessive structural distortions in response to the subsoil undergoing volume changes 
induced by a change in subsoil moisture. If foundations are not of a stiffness to resist or 
span the gaps caused by differential soil movements, and if the superstructure is not of a 
suitable flexibility, the foundation distortions resulting from the soil movements can be of 
a magnitude to cause excessive cracking of walls and distortion of internal fittings. This 
type of differential shrink and heave of soils give rise to two distinct distortion modes 
which cause damage to the superstructure. Fig 2 shows typical superstructure damage that 
would occur in the case of center lift. Center lift is a distortion mode resulting from the 
soil at the perimeter being more dry (perimeter soil may be shrinking) than the soil beneath 
the slab (soil beneath the slab is either at equilibrium, shrinking at a rate slower than the 
exterior soil, or heaving.) Fig 3 shows typical superstructure damage that would occur in 
the case of edge lift. Edge lift is a distortion mode resulting from the soil at the perimeter 
being more wet (perimeter soil may be heaving) than the soil beneath the slab (soil beneath 
the slab is either at equilibrium, heaving at a rate slower than the exterior soil, or 
shrinking.) The engineering profession lacks a universal method of analysis within a 
theoretical framework that utilizes fundamental properties of soil and footing behavior to 
enable a footing design that is satisfactorily accurate without resorting to empirical 
methods (Mitchell 1980). Since the late 1970’s strides have been made by the Post 
Tensioning Institute (PTI) to adopt a standard method for designing ground structures in 
expansive soils. The Post Tensioned Institute has adopted the methodology discussed 
herein and stipulates that post tensioned slab designs shall provide necessary rigidity to 
avoid damage to the buildings or houses when the supporting soil expands or contracts 
(PTI 1 996). However, a theoretical approach to predict vertical differential soil movement, 
as presented herein, has not been widely adopted by today’s practicing engineers. 
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Fig. 2. (Mitchell 1980) Center Lift Distortion Mode Fig. 3. (Mitchell 1980) Edge Lift Distortion Mode 

Expansive Soil, on a Molecular Level 

The Building Blocks of Expansive Soil 

Why are “expansive soils” so expansive? The answer can be found in the material 
on a molecular level. The magnitude of the surface area per mass and the molecular 
composition of expansive soils plays a major role in the soil’s ability to shrink and swell to 
such an impressive degree. “Expansive soils” are typically the fine grained soils as 
classified in the unified classification system. Fifty percent or more of a soil, by weight, in 
a sample that passes the number 200 sieve is classified as a fine grained soil. Fine grained 
soils are further broken down into categories of silts, clays and organics having either a 
high or low plasticity. The ability to break down and classify a soil according to the unified 
soil classification system provides an engineer with a wealth of information about the 
behavioral characteristics and how a particular soil may effect design. An in depth 
explanation of soil characteristics and the unified soil classification system can be found in 
Lambe and Whitman (1969). 

The following clay structural unit explanation is idealized and simplified and are 
presented to provide a picture of the basic behavioral characteristics that give rise to clay’s 
ability to undergo large changes in volume. The structural units that make up a clay 
particle are silicate sheets. Primarily, clays are made of sandwiches of: silica sheets, 
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gibbsite sheets, and brucite sheets. These sheets are combined in two distinct groups: 
two-layer sheet minerals, and three layer sheet minerals. The most important and most 
common two-layer sheet mineral is the kaolinite clay particle, see Fig. 4. Kaolinite is a 
two-layer mineral which consists of a silicate sheet bonded to a gibbsite sheet by 
hydrogen bonding and secondary valence forces. Another example of a two layer sheet 
mineral is serpentine which is formed by brucite bonding to silica, see Fig. 5. 
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Fig. 4. (Lambe and Whitman 1969) The Structure of Kaolinite 
(a) Atomic Structure (b) Symbolic Structure 



In the actual formation of these silicate sheets there may be some isomorphous 
substitution. Isomorphous substitution is when one kind of atom is substituted for 
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Fig. 5. (Lambe and Whitman 1969) The Structure of Serpentine 
(a) Atomic Structure (b) Symbolic Structure 

another atom that is very similar in structure and size. In other words aluminum may be 
substituted for silicon if there is a shortage of silicon and an abundance of aluminum 
during the formation stage of the clay mineral. This sort of substitution may give the 
mineral a net charge deficiency and may cause slight distortion in the crystal lattice due 
to the substituted atom’s size difference. Isomorphous substitution that occurs giving 
the clay mineral a large negative charge will give the mineral the capacity to attract 
polar charged water molecules. A double layer of molecules, such as water molecules, 
will try to form around a clay mineral to neutralize the negatively charged clay particle. 
This great affinity for water due to isomorphous substitutions provides an explanation 
as to why clay minerals, like montmoril Ionite, can have large volume changes at high 
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volumetric water contents. In Fig. 6 the clay particles are shown with the fully 
developed double layers they would have in pure water. 





Fig. 6. (Lambe and Whitman 1969) Soil Particles With Water and Ions 
(a) Sodium Montmorillonite (b) Sodium Kaolinite 

In kaolinite there is a very small amount of isomorphous substitution. However 
in the three layer sheet clay minerals there is much more isomorphous substitution. 
Three layer sheets are formed by placing a silica on the top of, and on the bottom of, 
either a gibbsite or brucite sheet. Fig. 7 shows the mineral pyrophyllite made of a 
gibbsite sheet sandwiched between two silica sheets. Fig. 8 shows the structure of the 
mineral muscovite, which is similar to pyrophyllite except that there has been 
isomorphous substitution of aluminum for silicon in muscovite. The net charge created 
by this substitution is balanced by potassium ions which serve to link the three layer 
sandwiches together. The most common three-layer minerals in soil are 
montmorillonite and illite. Montmorillonite is similar to pyrophyllite with the 
exception that there has been isomorphous substitution of magnesium for aluminum in 
the gibbsite sheets. Due to the isomorphous substitution the illite mineral has a 
moderately high capacity for volume change and montmorillonite has a very high 
capacity for volume change (Lambe and Whitman 1969). 
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Fig. 7. (Lambe and Whitman 1969) The Structure of Pyrophyllite 
(a) Atomic Structure (b) Symbolic Structure 
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Fig. 8. (Lambe and Whitman 1969) The Structure of Muscovite 
(a) Atomic Structure (b) Symbolic Structure 
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Gilgai 

Clay soil left to interact with surface climate conditions will develop an irregular 
surface microrelief pattern called gilgai, which in colloquial terms is referred to as “hog 
wallows.” This phenomenon is recognizable in areas where highly expansive soils are 
known to exist undisturbed for extended periods of time. These gilgai may be 
experienced while taking a very undulating ride: up the fairway of an old golf course, 
during farming evolutions in a native hay field, or down an old country road. As 
illustrated in Fig. 9 from Spotts (1974), the pattern develops in seven stages. First, the 
sedimentary clay is deposited. Following deposition, drying occurs, producing surface 
shrinkage cracks. Further desiccation results in the formation of major cracks 
penetrating far into the clay. The extent of penetration and spacing are determined by a 
number of factors of which the degree and rate of drying and the clay type are the most 
significant. Subsequent development of microrelief occurs as the clay is wetted and 
dried by cyclic variation of the climate. Because water penetrates deeper at the major 
cracks, greater swell occurs. As the next cycle of drying proceeds, the higher areas dry 
first, causing the major cracks to reappear in approximately the same locations. Thus, 
over a period of several years, an undulating surface pattern called gilgai develops on 
the soil surface (McKeen 1981). This phenomenon is important to realize since it is a 
natural occurrence of differential soil movement that must be understood by engineers 
so that appropriate models of this type of soil behavior can be properly applied to 
design. There are important parameters to consider when addressing volume strain that 
causes the gilgai phenomenon, such as: (1) The depth at which soil interacts with the 
climate, known as the depth to constant suction, z m . (2) The climate, or the annual 
wetting and drying cycles at the surface. (3) The type of surface conditions, such as; 
vegetation growth, bare soil, slope, or barriers of some type. (4) The soil characteristics, 
namely: suction compression index, y h (Zj), unsaturated permeability, k o (z s )|h c | , and soil 
diffusivity, a(Zj). The roles that each of these parameters play in predicting differential 
soil movement are presented. 
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Fig. 9. (Spotts 1974) Gilgai Stages 
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Suction, Matrix and Osmotic 

Soil suction is the measure of a soil’s affinity for water. Soil suction is an 
energy: water always moves down an energy gradient and it takes energy to cause 
volume strain in any soil. Soil scientists and arboriculture scientists refer to soil suction 
in terms of water potential. As defined by the 1960 Moisture Equilibria Symposium, 
soil suction is the negative gage pressure relative to the external gas pressure on the soil 
water to which a pool of pure water must be subjected in order to be in equilibrium 
through a semi-permeable membrane with the soil water. Suction is the term used 
principally by engineers for the thermodynamic quantity, Gibb’s free energy which is 
inherently negative, as seen in (1) below, and generates tension in the pore water 
stretching between soil particles (Lytton 1994). Soil suction is equal to 0.0 when the 



relative humidity is equal to 100%. A relative humidity value less than 100% in a soil 
would indicate the presence of suction in the soil. 

Components of Soil Suction 

Soil suction is quantified in terms of the relative humidity (1) and is commonly 
called “total suction.” “Total suction,” vp , has two components, namely: matrix suction 
(p), which is due to the attraction of water to the soil particle surfaces, and osmotic 
suction (7r), which is due to dissolved salts or other solutes in the pore water (2). 



Matrix suction will have the most effect on the volumetric changes in soil due to the fact 
that this component varies with the changes in the moisture environment. Osmotic 
suction will not change unless there is a change in the salt concentration within a given 
soil layer. 

Typical Suction Levels in the Field 

Fig. 10 illustrates the suction-vs-volumetric water content curve for a natural soil under 




( 1 ) 



V|/ = p + 71. 



( 2 ) 
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wetting and drying conditions. A common measure of suction is the pF-scale, in which 
pF is defined as (3), where p refers to the logarithmic value analagous to pH, and the F 
refers to free energy. 

Suction in pF = log 10 | h | (3) 

Where | h | equals the magnitude of suction in centimeters of water, a positive value. 
The various suction levels corresponding to field cases are listed below: 




Fig. 10. (Lytton 1994) Suction vs Volumetric Water Content Curve 

Field capacity (h^ = 2.0 pF) 

Wet Limit for clays (1^,= 2.5 pF) 

Plastic limit (3.5 pF) 

Wilting point of plants (h diy = 4.5 pF for vegetation) 

Tensile strength of confined water (5.3 pF) 

Air dry at 50% relative humidity (h dry = 6.0 pF for bare soil) 

Oven dry (7.0 pF) 
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Field capacity is an agricultural engineering term used to define the maximum amount 
of water that a soil can hold in the field. Fig. 10 shows field capacity to be a value of 
2.0 pF, but actual field values never quite drop below a value of 2.5 pF for clays (Lytton 
1995 and Mitchell 1980). Soil suction is a term that measures energy potential and is 
directly related to the soil’s potential for volume change. This information is referred to 
repeatedly throughout this document. 
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THE THEORY OF DIFFERENTIAL SOIL MOVEMENT 

Some notable breakthroughs have been made pertaining to the topic of 
predicting differential soil movement and are described in this chapter. The real 
challenge in geotechnical engineering, especially in the area of expansive soils, is to 
develop a consistent body of theoretical mechanics which adequately explains the 
behavior of the soil mass and to develop methods of measuring the relevant material 
properties critical for making appropriate design decisions. 

In the 1970’s the first attempt was made to design slabs-on-ground using a 
scientifically based theoretical approach. In 1973, Lytton and Woodbum prepared a 
report which chronicles a successful school building foundation design. The school 
foundation design procedure considered differentially expanding and shrinking clay due 
to the changing suction of the soil beneath the foundation. The suction profiles of the 
soil were predicted using a rational method which relied upon suction profiles measured 
on samples taken on the site at different times in the year. These suction values were 
used to predict the differential soil movement, and the foundations were designed based 
on the distortion modes caused by the magnitude and location of the expansive soil 
movements. Lytton noted that this differentially expanding and shrinking effect caused 
“mounds” to form beneath the slab. The foundation interacts with these mounds 
pressing down on the high spots and bridging the low spots. The foundation design 
procedures considers the limit cases in which the mounds produce two types of 
distortion patterns: edge lift pattern and center lift pattern, as shown in Fig. 1 1. These 
two distortion modes were used to design the slabs so that the loads of the 
superstructure would be carried across a span generated by edge lift, or carried in a 
cantilevered fashion out to the edge of the foundation in the case of center lift. 
Woodbum recorded the performance of the foundation systems after the designs were 
built and the school was put into service. Lytton and Woodbum reported post 
construction edge shrink values to be 2.75 inches which compared favorably with the 
predicted value of 3.1 inches. Just as Lytton predicted, the worst case differential 
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movement occurred in the first dry season subsequent to the construction of the slab 
during the wettest time of the year. In this post construction case the soil was “capped” 
and the moisture content was preserved in a 





Edge Heave 



Center Heave 



Fig. 11. (Lytton 1973) Edge Heave and Center Heave Distortion Modes 

state of heave under the slab while the perimeter soil shrank as it dried. This post 
construction case produced the worst case center lift mode. Subsequent weather cycles 
did not produce such extreme cases of differential soil movements as compared to the 
post construction case. However, subsequent weather cycles did produce differential 
soil movements as predicted. The model to predict this type of post construction 
differential soil movement is presented in the Post Construction section of this chapter. 

This design and research venture documented by Lytton and Woodbum’s 1973 
paper marked the beginning of a slab-on-ground design method that has been adopted 
by the Post Tensioning Institute as the standard for design of post-tensioned slabs-on- 
ground on expansive soils (PTI 1980, 1996). However, the objective of this research is 
to consider the methodology behind predicting the differential soil movement used to 
design such slabs. 

Soil Volume Strain 

Movements in expansive soils are generated by changes of suction which are 
brought about by the entry or loss of moisture or change in concentration of dissolved 
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salts in the pore water. The volume change that accompanies the change of suction (and 
water content) depends upon the total stress states that surround the soil. Within a soil 
mass, a decrease in the magnitude of suction results in an increase of water content. 

The volume of the soil also increases unless the surrounding pressure is sufficient to 
restrain the swelling. A conceptual graph of suction-vs-volume can be drawn using the 
relations of each to water content. This is illustrated in Fig. 12 on the plane 
corresponding to zero pressure. A similar graph can be drawn relating pressure (total 
stress) - versus - volume on the plane corresponding to zero suction. The simultaneous 
change of the magnitude of suction (decrease) and pressure (increase) results in a small 
change of volume, following the path from Point A to Point C on the pressure - suction 
- volume surface. The magnitude of suction decreases from Point A to Point B while 
the pressure increases from Point B to Point C. The volume change process can be 
viewed as the net result of two processes; 

a. Increase of volume from A to B at constant mechanical pressure or total 
stress. 

b. Decrease of volume from B to C at constant suction. 



Volume 




Fig. 12. (Lytton 1994) Mean Principal Stress-Volume-Matrix Suction Diagram 
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For small increments of volume change on this surface, the volume strain, is linearly 
related to the logarithms of both pressure and |suction| . The soil volume strain equation 
which shows the relation between matrix suction, mean principal stress, and osmotic 
suction is (4), (Lytton 1977, 1994). 
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In the procedures described herein the assumption will be made that the osmotic 
component has no change. If the case exists that the osmotic component has no real 
change, which is the case most of the time, then that term within the volume strain 
equation can be neglected. The suction compression index , y h (Zi), and the mean 
principal stress compression index, y 0 (z i ), will be discussed in more detail in the section 
titled “The Suction/Mean Principal Stress Compression Indexes”. The initial and final 
suction profiles, h,(Zj) and h F (Zj), will be discussed in more detail in the section titled 
“Suction Profiles”. The initial mean principal stress, ct„ is the stress level below which 
there is no suppression of volume strain by overburden pressure. This is the stress at 
which the curve in the pressure-volume plane in Fig 12 begins to depart from horizontal 
as pressure increases from zero and has been observed to be approximately 40 
centimeters of soil (Lytton 1994). The final stress at depth Zj a F (Zj), is the mean 
principal stress at a depth z, at the time the volume strain is desired and is given by (5). 
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( 5 ) 



where k„ is: 



k 0 = 0.00 when the soil is badly cracked. 
k 0 = 0.33 when the soil is drying. 
k 0 = 0.67 when the soil is wetting. 

k 0 = 1 .00 when the cracks are closed and the soil is swelling. 
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Once the volume strain for an incremental layer of soil is known the vertical strain for 
that layer is estimated using a cracked fabric factor, f ( , as shown in (6). The cracked 
fabric factor is the percent soil volume strain that is directed vertically in a soil column. 
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v v / 



McKeen (1981) back calculated vales of f,and found that f ( = 0.5 when soils are drying 
and fj = 0.8 when soils are wetting. The total vertical differential movement in a soil 
mass, Az , is then found by summing the product of vertical strain and the incremental 
depth to which they apply as shown in (7). 



Az = X 




• increment 



( 7 ) 



The volume strain theory previously described was developed by Lytton (1973, 1977, 
1994, and 1995). All symbols are defined in the List of Symbols. 



A Simplified Method for Identifying Predominant Clays 

There are three areas of soil characteristics that must be known before the 
volume strain theory can be applied, namely: the soil suction compression index, the 
soil suction profiles, and the mean principal stresses within the soil. The mean principal 
stress was described earlier and suction profiles will be discussed in a later section. 
However, the soil suction compression index, as described in the next section, depends 
on proper identification of the predominant clays in soils. A simplified method for 
identifying these predominant clays is discussed in this section. 

Before Lytton developed the methodology to predict vertical differential soil 
movement, Pearring (1968) and Holt (1969) had completed a correlation chart to aid in 
the identification of predominant clay mineral in a given soil. According to Pearring 
and Holt’s research, the two parameters used to classify the clay minerals are the cation 
exchange activity, CEAc(Zj), and activity ratio, Ac(Zj). Obtaining the measure of these 
two parameters requires the plasticity index, PI(Zj), the cation exchange capacity, 
CEC(Zj), and the percent of clay in the soil passing the number 200 sieve, Clay(Zj), (8), 
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(9), (10), (11). 



PIfe) = LL(z i ) - PL(Zj) 



( 8 ) 




(9) 




( 10 ) 




( 11 ) 



All symbols are defined in the List of Symbols. Once the cation exchange activity and 
the activity ratio are known for a given clay mineral the type of clay is determined fron 
Fig 13. 



The difficulty comes in arriving at the measure of cation exchange capacity. 
Until Pearring and Holt’s method of identifying predominant clays came to light the 
geotechnical engineer had to rely on experience or the more reliable, expensive, and 
arcane techniques used by clay mineralogist and soil physicists. However, Pearring 
and Holt’s method still required engineers to find the cation exchange capacity (CEC) 
of the soil, which is a soil property not typically evaluated in normal soils test 
laboratories. Mojekwu (1979) realized this problem and developed a simplified way to 
arrive at the CEC. His work involved data collected for a wide range and number of 
soil samples taken in the state of Texas. Among other data, Mojekwu collected and 
analyzed the Atterberg limits of these soils and found that very good correlations could 
be made that linked the plastic limit and the liquid limit to the CEC of a specific soil. 
The method he developed requires the simple Atterberg limits tests to find the plastic 
limit, PL, and the liquid limit, LL. Mojekwu ran a simple regression analysis, the 
dependent variable being CEC, the independent variables being PL and LL respectively, 
and produced very high R squared values of .9941 and .9942 for (12) and (13) 
respectively. 



CEC(Zj) - PL(Zj) 117 
CEC(Zj) = LL(Zj) 0 - 912 



( 12 ) 

( 13 ) 
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The larger of these two approximated values of CEC is used to find the CEAc, which is 
then used along with Ac to determine the predominant clay minerals in Fig. 13. These 
findings by Pearring, Holt, and Mojekwu paved the way for research by McKeen in 
developing a simplified method to identify the suction compression index, which is 
critical to predicting differential soil movements. 
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Fig. 13. (Pearring 1968 and Holt 1969) Chart to Identify Predominant Clays 



The Suction/Mean Principal Stress Compression Indexes 

Suction forces on a soil element are body forces rather than externally applied 
forces at the surface of the element. Because the suction stresses are dominating the 
effective stresses in expansive soils, volume change studies are made using total suction 
values. It is the slope of the volume change versus suction curve that quantifies soil 
response to moisture changes. Thus, the suction compression index, y h (Zj), is defined as 
the slope of the volume-total suction curve, see Fig 14 (Lytton 1977). McKeen (1981) 
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analyzed soils data in which he measured volume strain corresponding to soil suction 
values to arrive at a corresponding suction compression index for particular types of 
expansive soils. In McKeen’s study, he measured the soil volume change using the 
conventional oedometer for cases when the soil was un-cracked and using the 
coefficient of linear extensibility (COLE) method when the clay was cracked. The 
measurements of soil suction were taken utilizing a variety of methods as well. The 
primary methods were the thermocouple psychrometer and filter paper methods. Using 
the basic concept 



described above, McKeen developed the chart method to determine y h by analyzing a 
large number of Soil Survey samples, adopting the methodology for determining 
predominant clays by Pearring (1968), Holt (1969) and Mojekwu (1979), and produced 
Fig 1 5 which gives values of y hl00 without requiring suction tests. The value of y hl00 
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Fig. 14. (Lytton 1977) Pressure-Volume-Suction Surface 
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typically ranges from .033 to .220. This method requires only the activity ratio and the 
cation exchange activity of the soil to find the corresponding y h for particular soils. The 
chart method is used in the procedures to calculate y h in the software design format 
procedures. 

Another factor to consider in calculating soil volume strain is the effect of 
applied loads in suppressing the suction volume strain by overburden pressure. The 
mean principal stress compression index is required to complete the term to calculate 
the amount of suppression to the suction volume strain by overburden pressure. The 
mean principal stress compression index is calculated by (14) (Lytton 1977). 

V V 
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Fig. 15. (McKeen 1981) Chart to Determine the Suction Compression Index 



The mean principal compression index is also calculated by the commonly used 
compression index, C c , by (15) (Lytton 1994). 




o 



( 15 ) 



Typical values of y a were calculated from results of swell tests and swell pressure tests 
reported by Vigayvergiya and Ghazzaly (1973) according to (14). The range of these 
calculated values of y CT fell between 0.02 and 0.09 indicating that the ratio of y h /y a is 
typically at or slightly below 1.0. Therefore it is slightly conservative but convenient to 
assume y h is equal to y CT Until there is research that develops a more practical and cost 
effective method for obtaining more accurate values of y a , that coefficient will assume 
the value of y h for calculations of volume strain in the software design format 
procedures described in this document (Vigayvergiya and Ghazzaly 1973). 

Estimates of Unsaturated Soil Properties 

Before a clear explanation of the remaining sections is possible, estimates of 
unsaturated soil properties must be defined. Specifically, definitions are required for the 
following soil properties, namely: the slope of the suction versus volumetric water 
content wetting line shown in Fig. 10, S(Zj), Mitchell’s diffusion coefficient, a(Zj), and 
Mitchell’s unsaturated permeability, k o (z ; )|h 0 | orP(Z|). The slope of the suction versus 
volumetric water content wetting line is given by (16). Mitchell’s diffusion coefficient 
is given by (17). Finally, Mitchell’s unsaturated permeability is given by (18). The 
process flow sheet that describes procedures for calculating all soil properties required 
for the software design can be found in Appendix A-88. 

S(Zj) = -20.29 + 0.1 555 ■ LUz; ) - 0.1 1 7 • PI( Zi ) + 0.0684 • v( Zj ) (16) 

a(Zj) = 0.0029 - 0.000 1 62 • S(Zj ) - 0.0 1 22 • y „ (z t ) (17) 



The above equations were developed by Jayatilaka et al. (1992) and Lytton (1994). All 
symbols are defined in the list of symbols. 




(18) 
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Edge Moisture Variation for Center Lift and Edge Lift Cases 

As discussed earlier, the effects of vertical differential soil movement on 
commercial and residential types of slabs-on-ground is dictated by the perimeter 
moisture effect cases and the total amount of vertical movement between two distinct 
locations separated by a characteristic distance. The characteristic distance is described 
as the edge moisture variation distance. At some time after construction, the soil under 
the slab will reach an equilibrium moisture content. However, the soil at the edge will 
vary in moisture content dependent on the local annual weather cycle. Additionally, the 
soil underneath the slab will vary in moisture content up to a distinct distance, the edge 
moisture variation distance, that is controlled by the permeability and diffusivity of the 
soil. This differential in moisture content of soil underneath the slab and the edge of the 
slab is what causes the center lift and edge lift distortion modes. As part of research 
conducted by Jayatilaka et al. (1992) a calibrated finite element program with coupled 
transient moisture flow and elasticity that had been used in the study of vertical 
moisture barriers provided an ideal means to study the edge moisture variation distance. 
The unsaturated soil properties discussed in the previous section in addition to the 
Thomthwaite Moisture Index were used to determine the relation of the edge moisture 
variation distance. 

The Thomthwaite Moisture Index (TMI), which measures the location’s 
potential for evapotranspiration, is defined as a climatic index that measures the water 
balance of an area to determine if there is a water deficiency or water excess. 

Both edge lift and center lift conditions were explored using several hundred 
runs with the program (Jayatilaka et al. 1992). Center lift conditions were simulated by 
a one year dry spell following a wet suction profile condition. Edge lift conditions 
were simulated by a one year wet spell following a dry suction profile condition. The 
edge moisture variation distance was considered to be that distance between the edge of 
the foundation and the point beneath the covered area where the suction changed no 
more that .2 pF during the entire period of simulation (Lytton 1994). 
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The dry and wet conditions used annual suction variation patterns that were 
appropriate for each of nine different climatic zones ranging from a Thornthwaite 
Moisture Index of -46.5 to +26.8, spanning the range found in Texas. The resulting 
edge moisture variation distances are shown in Figs. 16 and 17. No distance less that 
2.0 feet should be considered for design purposes (Lytton 1994). In the procedures for 
the software design, the edge moisture variation distances, e me (Zj) and e mc (Zj), are 
calculated using Figs. 16 and 17. 

In 1 997 Lytton developed a method to determine the edge moisture variation 
distance for a particular design return period e mr . In designing slabs-on-ground 

recognition must be taken of the length of the time these structures must be in service, 
and of the severity of the weather patterns that may occur during the expected life of the 
structure. The return period of hydrologic events is appropriate to use in estimating the 
design criteria for slabs-on-ground on unsaturated soils. Determining the edge moisture 
variation distance for a selected return period involves applying an equation that relates 
the edge moisture variation values for the 10-year and 50-year return periods to the 
design edge moisture variation value for the selected return period, (19). 



e„ = e m + e m -e 



( e m 5 „ e m„, )' 



Z r Z.O 



k Z 50 Z I0' 



(19) 



The edge moisture variation distance for a 10 year return period can be obtained from 
Fig. 1 8. These edge moisture variation distances were derived by back calculation from 
slabs which were performing successfully in San Antonio, Dallas, and Houston. None 
of the slabs were more than 10 years old at the time. It can be argued that the design 
values of edge moisture variation distance represent a 10 year return period, e . The 
edge moisture variation distance charts shown in Fig. 16 and Fig. 17 represent a 50 year 
return period, e . The use of the Gumbel probability density function, which is 
commonly used to represent the probability of weather events, may be used to establish 
the risk level that is desired for design in accordance with the expected service life. 

The z, scores used in (19) are computed from the Gumbel cumulative probability 
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distribution curve, (20). The return period is represented by the variable r. p and p are 



z, = 



( 20 ) 



- in "-L 



both shape factors and can be assumed as one. Solving for an edge moisture variation 
distance for a particular return period between 10 and 50 years is a simple two step 
process: (1) simply substitute the return periods of 10, 50 and the selected return 
period r into (20) and solve for z 10 , z 50 and z, respectively, and (2) Substitute values of 
e from Fig. 18, e from Fig. 16 or Fig. 17, and values of z 10 , z 50 , and z, previously 



calculated into equation (19) and solve for the e, n . A common design period for 
residential and light commercial structures is r = 20 years (5 percent risk). 
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Fig. 16. (Lytton 1994) Edge Moisture Variation Distance for Center Lift Conditions (50 yrs) 
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Fig. 18. (Wray 1978) Edge Moisture Variation Distances for Edge/Center Lift Conditions (lOyrs) 
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Equilibrium Suction 

Before the equilibrium suction profiles, h^Zj), wet and dry limit suction profiles, 
h mjn (z i ) and h^fo), transient state suction profiles, h(z,,t), and depth to constant suction, 
z, n can be generated, the equilibrium suction , h in , must be known for particular soil 
profiles and surface conditions. Gay (1994) did extensive work in the area of 
calculating the mean volumetric moisture content for a given soil mass dependent upon 
the soil’s depth of available moisture, d am , and the location’s potential 
evapotranspiration. The depth of available moisture, d am , is defined as the maximum 
depth of moisture available for use by transpiring vegetation, which is stored within the 
soil zone down to the depth to constant suction. The Thomthwaite Moisture Index is 
used to ascertain if a climate is arid or humid (Thomthwaite 1948). Gay’s theory 
considers d am and TMI. With these, he developed a set of functional relationships that 
are used to calculate the mean volumetric moisture content, which are then used to 
calculate the equilibrium suction value for a particular soil profile and location. The 
approach he used to solve for these equations was developed by Juarez-Badillo (1975) 
in which functional domains for the problem are first established as shown in Fig. 19. 
Function (21) was developed to satisfy the boundary conditions in Fig. 19. Then 
through an assumption of linear proportionality of the rates of change of the two 
complete functions, Gay obtained, (22) and (23). 



f(d) = 





( 21 ) 



dT df(d) 
Y T " f(d) 




= In 



1 




d^ -d 


^ am 


1 


1 


^ am _ ^ 1 


d a«J 



( 22 ) 



(23) 
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Fig. 19. (Gay 1994) Domains of Functions f(d) and g(T) 



Solving (23) for d gives (24). 



d = 




(24) 



All symbols are defined in the list of symbols. Gay then determined how the 
coefficients y, d,, and T, depend upon d am using the pattern search technique and arrived 
at the expressions (25), (26), and (27). 



y = 0.039337d m + 1.357033 


(25) 


d, = 0.449079d am + 0.304560 


(26) 


T, = 0.06265 ld am + 59.53593 


(27) 



These relationships facilitate the calculation of mean moisture depths, d m , for all values 
of TMI for any value of available moisture depth between 10cm and 50cm, as illustrated 
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in Fig. 20 (Gay 1994). The development of (24) is important because it is used to solve 
for the mean volumetric water content, 0 m The value 0 m is applied as a target to solve 
for h m in an iterative process using Nieber’s (1981) equation (34). Nieber’s equation 
relates suction to volumetric water content. The method to apply this iterative process 
to solve for the equilibrium suction value can be found in Appendix A-90. 




50 cm 
40 cm 
30 cm 
20 cm 
10 cm 



Fig. 20. (Gay 1994) Theoretical Relationships for Mean Annual Moisture Depth vs TMI 



The equation used to solve for the mean annual moisture depth, (24), developed 
by Gay is rearranged to arrive at the mean volumetric water content, (32). Below is the 
step by step process to obtain the equation to solve for 0 m . First note the following 
equations from Gay (1994), (28), (29), (30), (31). Fig 21 provides a graphical 
explanation for some of the symbols in these equations. All symbols are defined in the 
list of symbols. 

d am 

Z ° 0 w! -0*jr 



( 28 ) 
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e... = 



d m +d d 0 ' 



A _ d wl 

e„ - 



®diy _ 



l dry 



Dividing through equation (24) by z c gives (32). 



d. 

z„ 



1 + 



d -d, 

am 1 



yV 



vT,/ j 

Solving for dn/z,. in (29) and substituting in (32) gives (33). 



(29) 

(30) 

(31) 



(32) 



8.- / 9wl ~ 9d,? n +9 



dry 



(33) 



1 I d am d l 

T y 
d| y 

The steps to arrive at equations which solve for 0 wl , 0 dry , and 0 m using Nieber’s 
equation (34) are given below. 







fe s -0 r V 
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CD 



(34) 



Introduce dummy variable r, (35). 

kf ^-0 

r A 0 - 0 r 



( 35 ) 
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Solving (34) for 0 gives (36). 



0 = 



Qs + r-Qr 
1 + r 



(36) 



Substituting 





or 




back into (36) for r gives (37), (38), or (39) 



depending on which volumetric moisture content is to be found. 



MOISTURE CHANGES IN UNSATURATED SOILS 



pF 5.7 pF 4.2 pF 1. 6-3.2 pF 1.0-1. 5 




df c = depth of moisture at field capacity 
dam = depth of available moisture 
Z r = rooting depth in soil 



Fig. 21. (Lytton 1994) Moisture Changes in Unsaturated Soils Between Typical Wet and Dry States 
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(37) 



(38) 



(39) 



The equilibrium suction value, can then be solved through an iterative process using 
equations (33), (37), (38) and (39) as described in Appendix A-90. The equilibrium 
suction profile, h, n (Zj) for any pervious soil profile is defined by decreasing the 
equilibrium suction value at the surface, h m a centimeter in suction for every centimeter 
increase in depth (McQueen and Miller 1968). The process for determining the 
equilibrium suction profile is detailed in Appendix A-96. 



Suction Profiles 

It is important to note that moisture changes within a natural soil are determined 
by the ratio of the period of evaporation to the period of precipitation. A measure of 
this period of evaporation to the period of precipitation is determined through a method 
developed by Thomthwaite (1948). In the previous section the equilibrium suction is 
determined based on relations to the climatic water balance expressed in terms of the 
Thomthwaite Moisture Index (1948). When a marked separation occurs between wet 
and dry seasons, a large seasonal variation in soil moisture content occurs, whereas in 
areas which are either predominantly dry or predominantly wet for most part of the year, 
the changes in soil moisture content is not so marked. A location where predominant 
clays are abundant in which the climate is marked by this large separation in the periods 
of evaporation and periods of precipitation will experience problems with differential 
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movement to the greatest degree. The state of moisture is dependent on the state of 
suction and tends towards an equilibrium suction profile, as defined in the previous 
section, with a moisture source either at the boundaries or within the soil mass. 

When the state of suction is measured at intervals of depth in a soil profile, the resulting 
relation between suction and depth is termed the suction profile. Mitchell measured 
actual suction profiles in the area of Adelaide, South Australia and found a very 
consistent and typical trend that could be modeled. In his study he developed a method 
based on constitutive soil properties to model the suction profile within a soil profile 
base on the local annual weather cycle, n, soil diffusion, ct(Zj), equilibrium suction 
value, hn,, and the amplitude suction value, U G The maximum limit suction profile, 
h max (z i ), and minimum limit suction profile, h^z,) begin at the limiting field condition 
suction value at the surface, and decrease at an exponential rate toward the equilibrium 
profile, h^Zj). These limit suction profiles will approach the equilibrium suction profile 
at the depth to constant suction, z, n , which will be defined in the next section. The limit 
suction value depends on the type of surface conditions and is defined in the section 
titled “Typical Suction Levels in the Field”. The dry limit suction value is usually 4.5 
pF for soil with vegetation at the surface and 6.0 pF for bare soil at the surface. Refer to 
Fig. 22 for a typical maximum limit suction profile. The wet limit suction value, which 
is somewhat drier than the field capacity suction value, is usually 2.5 pF. Refer to Fig. 
23 for a typical minimum limit suction profile. Mitchell uses an idealized model of the 
annual weather surface cycle which he represents by (40). This is a good representation 
for design because it assumes that soil at the surface approaches the maximum and 
minimum moisture conditions at two distinct and separate times of the year. As stated 
earlier, this separation between wet and dry periods produces worst case differential soil 
movement scenarios. Mitchell combines the soil suction relationship with the 
relationship that describes the local surface weather cycle and arrives at equation (42). 

Mitchell found that the suction change due to the effects of climate, drainage and 
site cover is a periodic function of time. The suction at any time and at any depth in the 
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soil profile is determined by solving the diffusion equation, (41), for a particular 
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Fig. 22. Typical Dry Limit Suction Profile for Bare Soil at the Surface 




Fig. 23. Typical Wet Limit Suction Profile for Bare Soil at the Surface 
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set of boundary conditions. The solution to arrive at the soil suction profile equation 
(42) can be found in Mitchell’s 1981 report (values of h(Zj,t), h^z^t), and U 0 are all in 
units of pF). 



h(Zj = 0,t) = cos 



27tn-t-(Zj =0) 



n-7i 



a(Zj = 0), 



(40) 



dh d 2 h 
8t~ a dz 2 



h(z i5 t) = h m (z i ) + U 0 e 



t * I n-7i 



COS 



27t •n*t-(z i )- 



n*7t 

a(Zi). 



(41) 

(42) 



Based on the surface and depth boundary conditions, referred to as moisture effect 
cases, (42) is applied to produce suction profiles for many surface conditions and soil 
profiles. Appendix A contains detailed algorithms that describe precisely how to apply 
(42) to arrive at suction profiles resulting from combinations of different moisture 
effect cases and soil profiles. Later sections in “The Procedures, in a Software Design 
Format” chapter will provide a more detailed discussion with respect to application of 
(42) in producing suction profiles. 

Depth to Constant Suction 

The depth to constant suction is the depth in a soil profile to which there is no 
longer a significant seasonal change in soil suction. This depth has been reported to be 
anywhere from 6 feet to as much as 26 feet dependent upon the type of soil, seasonal 
weather trends, and the type of surface conditions such as vegetation and slope. Mohan 
and Rao reported where the active zone extended to a depth of 13 to 16 feet, in the 
black cotton soils of India. Kassiff reported movement of an exposed clay at depths up 
to 20 feet in Israel, with the largest differential movements being measured within an 
active zone of 6 Vi feet, where moisture changes are maximal (Mitchell 1980). As one 
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uses the procedures contained herein to predict the differential movement in soils, it 
becomes clear that most of the movement does in fact occur nearer the surface. This 
reduction in the amount of differential soil movement with an increase in depth is 
attributable to the combination of small seasonal changes in soil moisture content with 
an increase in depth, and the increase in restriction to volume strain due to overburden 
stresses with an increase in depth. The algorithm to solve for the depth to constant 
suction is based primarily on Mitchell’s equation (42) that is described in the previous 
section. 

Lytton and the author are responsible for developing the following steps to 
arrive at the equation to calculate the depth to constant suction. Only the limit profile 
aspect of Mitchell’s equation, (42), will be used. Setting the cosine term equal to one 
forces this equation to generate only the maximum and or minimum suction envelopes. 
Refer to (43). 



volume strain for post equilibrium moisture effect cases. Equation (43) is rearranged to 
isolate the exponent term, (44). 




1.0 



h ( z i>t) = M z i) + U o e 



(43) 



Mitchell’s equation in the form of (43) is used in most of the algorithms to calculate the 



h( z ijt)- h m (Zj) = U 0 e 




(44) 



(44) is then divided through by U 0 , (45). 




U 



o 



( 45 ) 
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The last step is to take the natural log of both sides and solve for z,, (46). 



l u o 



/ 




(46) 



The depth to constant suction, z^, has been obtained when the amplitude at depth z ; is 
equal to 0.2 pF. Later sections in “The Procedures, in a Software Design Format” 
chapter will provide a more detailed discussion with respect to application of (45) and 
(46) in calculating the depth to constant suction. 

Horizontal Velocity in Expansive Soils 

The ability to calculate horizontal velocity of water flow in soils is critical for 
the application of the soil volume strain theory. Predicting vertical differential soil 
movement when a horizontal barrier is employed, or when a moisture effect case is 
located a distinct distance from the slab’s edge depends upon the soil’s horizontal water 
flow characteristics. The horizontal velocity is calculated and then used to calculate a 
suction profile a distinct distance from another known suction profile. Lytton and the 
author developed the horizontal velocity equation by applying Darcy’s law, (47), to the 
unsaturated permeability relationship for clay soils developed by Mitchell, (18), to 
arrive at the horizontal velocity equation for unsaturated soils, (52). The steps to derive 
(52) are described below. First Darcy’s law, (47), is solved for dh, (48). 




(47) 




(48) 



Darcy’s permeability constant is set equal to Mitchell’s unsaturated permeability, (49). 
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k = 



K 0 (z,).h 0 



(49) 



Substitute (49) in for Darcy’s permeability in (48) and separate variables dependent on 
distance and variables dependent on suction to arrive at (50). 



V,(z.) d = dh 

K 0 ( Zi )h 0 X h 



(50) 



Integrate the left side of (50) from zero to a distance x, and the right side of (50) from an 
entry suction value, h entry (Zj), to an exit suction value, h exit (Zj). Refer to equation (51). 

V x (z.) 



K 0 (z.) 



i dx= 1 

' } “o o h oltty ( Zi ) 



dh 

h 



(51) 



Solve equation (51) for V x (Zj) to arrive at equation (52). 



V x (z f ) = 



K 0 ( Zi )h 0 



•In 



/ h ^ 



‘entry (^i ) 
h exit (^i ) • 



(52) 



The entry suction value is the suction value at x = 0. The exit suction value is the 
suction value at x. Once the horizontal velocity is known between two distinct locations 
in an incremental soil layer, the suction value for any point between those two distinct 
locations is solved for using equation (53). Equation (53) is derived by solving for the 
exit suction in (52). 



V x ( Z j)x 



h .xi.( Z i)= h e„ W ( Z i)- e K ° (Zi) "° 



(53) 



Post Construction Theory 

The worst case differential soil movement may occur in the first dry season 
subsequent to the construction of the slab during the wettest time of the year, or 
during the first wet season subsequent to construction of the slab during the driest time 
of the year. This type of moisture effect case is referred to as the post construction 
case. The post construction case is a transient case for calculating volume strain of 



42 



expansive soils immediately after construction and before the soil under the center of 
the slab has reached an equilibrium moisture content. For example, when a slab is 
constructed the soil is “capped” and the moisture content changes slowly toward its 
equilibrium value while the perimeter soil shrinks and swells with the change in 
seasonal moisture. This type of scenario provides conditions for the post construction 
moisture effect case. The moisture content of the “capped” soil looses or gains 
moisture for some period of time until the soil under the center of the slab reaches an 
equilibrium moisture content. The rate at which the soil approaches equilibrium is 
dictated by the soil’s diffusivity and unsaturated permeability. Once the soil under the 
slab reaches the equilibrium moisture content, subsequent suction distributions and 
differential movements must be determined using the post equilibrium methodology. 

The “ Damped” Slab Suction Profile 

Lytton and the author developed the following set of equations to model the soil 
under the slab subsequent to construction in a transient case in which a soil suction 
profile approaches an equilibrium suction profile. The idea is to calculate the soil 
suction profile under the slab at the time of construction, h s)ab (Zj,t const/dry ), and then 
dampen this profile until it approaches the equilibrium profile. The soil suction profiles 
under the slab and at the slab edge are solved, with respect to time, based on some 
knowledge of the number of annual weather cycles. The annual weather cycle will be 
assumed to vary according to the function described by (40). To demonstrate this 
method, the annual weather cycle is assumed to be one, n = 1. This means that, in one 
year there will be one wet season and one dry season. The driest, or wettest, month of 
the year must be known. In this case, June is selected as the driest month of the year. 
Fig. 24 is an example of a weather cycle modeled using (40) for bare soil at the surface 
(n = 1, and June is the driest month of the year). 

The approach used to derive the equations to damp the suction profile at time of 
construction back to the equilibrium suction profile was developed by Juarez-Badillo 
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(1975) in which functional domains for the problem are first established as shown in 
Fig. 25. Function (54) was developed to satisfy the boundary conditions in Fig. 25. 

For Bare Soil At The Surface 




Fig. 24. A Weather Cycle Modeled by a Cosine Function, n = I, June is the Driest Month 



f (h) = - (Case 1 , when h,,, is greater than h const .) (54) 

1 m 1 in 1 const 

f(h) = 7 1 7 ~ 7 T (Case 2, when h, n is less than h const .) (54) 

^ ^ const 

The suction value with respect to time and depth for a profile that is modeled to tend 
toward an equilibrium suction was solved for cases 1 and 2. Both cases arrived at the 
same expression (59). The steps to arrive at (59) are described below. Through an 
assumption of linear proportionality of the rates of change of the function equations 
(55) is established, and the constant of proportionality, p, is equal to 1. 

df (h) _ dt 
f(h) P t 



( 55 ) 
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Fig. 25. Domains of Functions h, t, and f(h) 



The left side of equation (55) is integrated from an arbitrarily chosen suction value, h, to 
a suction value half way between h m and h const , h, The right side of equation (55) is 
integrated from a time corresponding to h„ t„ to the time required to approach h, t. This 
integration is given by equation (56). 



1 



1 



In h 1,1 ~ h h --W = p . ln l 

1 L_ t, 



(56) 



h m -h, h m - h const 

Equation (56) is rearranged to collect all constant variables on one side of the equation 
and set equal to a variable, r. Refer to equation (57). 
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h . ~ h n 
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= r 



(57) 



Now equation (57) is rearranged to solve for the suction value, h, given by (58). 



h = 



hcons, +r ‘ h n 

1 + r 



(58) 



Substitute (57) back in to (58) to arrive at the equation used to solve for the suction 
profile under the slab some time, t, after the slab was “capped” as a result of slab 
construction, (59). 
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h(z, t) = 




(59) 



The Time Equation 

Before equation (59) can be applied, a relationship must be derived to determine 
the amount of time required for water to move upward or downward to change the 
suction at any given depth half way between its initial value and the equilibrium value 
of suction. This time equation for t, is given by (61). The derivation of (61) is as 
follows. Darcy’s law of permeability is the relationship used to introduce the time 
variable. Velocity is defined as a distance through which water flows in soil, x, traveled 
in some amount of time, t,. Velocity described in terms of permeability of water in 
expansive soils in which the time, t,, required to obtain half of the change in suction (the 
time to obtain half of the change in suction was arbitrarily chosen), 0.5 • Ah , is given by 



Equation (60) is solved for t„ (61) In equation (60), the variable x is replaced with 
equation (44), where the ratio of suction values is set equal to one half. 



The time variable represented by (61) is used in (59) to solve for the suction profile 
some time, t, after soil is “capped” as a result of slab construction. 



(60) 




(60) 




t 




(61) 
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THE PROCEDURES, IN A SOFTWARE DESIGN FORMAT 

Introduction to the Software Design Approach 

The tools described in “The Theory of Differential Soil Movement” chapter are 
all that is needed to develop the equations and terms used in the procedures described in 
this chapter. Appendix B is a graphic representation suggested as a flow of vertical 
differential soil movement algorithms in a windows environment type software design 
layout. The suggested software design layout, Appendix B, is modeled after the 
PTISLAB WIN 1.0 software package owned by Geostructural Tool Kit, Inc. (1996). 
The focus of this document is not software design or marketing, those skills are left to 
experts in that area of work. The focus is to present the equations and terms in easy to 
follow algorithms, as described in Appendix A, so that engineers or programmers can 
apply this soil volume strain technology in everyday practice. 

A question that a practicing engineer may ask is, “What practical engineering 
cases can these procedures address?” The answer is, cases involving a surface drained 
soil profile in which the slope is stable and surface conditions are those described 
below: 

1 . Bare soil at the surface. 

2. Grass at the surface with shallow roots. 

3. A flower bed at the surface with a known depth of flower bed zone, d zone . 

4. A tree with a known depth of root zone, d zone . 

The procedures consider design effects that are introduced to change the magnitude of 
the vertical differential soil movement. These design effects include: vertical barriers, 
horizontal barriers, and the effects of locating the surface moisture effect a distinct 
distance from the foundation, d source . All of these design effects can be rigorously 
analyzed to enable an engineer the ability to choose the effect which will provide the 
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desired design results. 

The algorithms contained herein address all cases and design effects listed above 
by way of two design approaches: (1) the post construction approach, and (2) the post 
equilibrium approach. An engineer would apply the algorithms, used in both 
approaches, to address all possible combinations of cases and design effects presented. 
The engineer compares the results of these cases to find scenarios that produce vertical 
differential soil movement that generates maximum bending, shear and deflection. This 
information is used to generate adequate foundation designs. The most promising 
scenario is to use the power of today’s computing technology to perform the iterative 
calculations and tests required to arrive at these design quantities. A programmer could 
apply software technology, such as PTISLAB WIN 1.0, to model a slab’s interaction 
with the differential soil movements produced by the algorithms contained herein to 
arrive at these maximum design cases. 

The algorithms are presented in a format that can be easily programmed using 
simple spreadsheet type software. First, input required for all calculations are entered. 
Soil properties and other related soil characteristics are calculated and stored for use in 
algorithms presented later. Equilibrium suction values, h^ and depth to constant 
suction values, z^, are calculated for a specific location and specific sets of surface 
conditions. With the above complete, all moisture effect cases and design effects listed 
previously are addressed using the post construction solutions and the post equilibrium 
solutions. 

Input 

When considering all cases and all calculations required to complete algorithms, 
an accounting of the variables required for input becomes important. Great effort is 
made to keep the complex nature of data required for input down to a minimum. The 
goal, relative to input of data, is to ensure that all variables required for input are items 
that are easy and economical for practicing engineers to obtain. This is the same goal 
authors had in mind while developing the simplified method for identifying 
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predominant clay minerals (Pearring 1968, Holt 1969 and Mojekwu 1979), the methods 
for obtaining the suction compression index and the mean principal stress compression 
index (McKeen 1981 and Viayvergiya and Ghazzaly 1973), and the method for 
obtaining the edge moisture variation distance (Jayatilaka et al. 1992). Below is a full 
list of items required for input to perform every calculation in the algorithms described 
in Appendix A. 

A,B Gardner’s coefficients (Appendix A-90 and Table 2) 

B Number of boring samples (Appendix A-84) 

D c Depth to the c th suction measurement in a sample (Appendix A-87) 

L s Depth to the bottom of layer s (Appendix A-84) 

LL S Liquid limit for layer s, in percent (Appendix A-84) 

M b Number of layers in the b 0 ’ boring sample (Appendix A-84) 

N Number of suction measurements in a sample (Appendix A-87) 

PL S Plastic limit for layer s, in percent (Appendix A-84) 

TMI Local annual Thomthwaite Moisture Index, (Fig. 26 and Appendix A-83) 

X b X coordinate for the b‘ h boring sample (Appendix A-84) 

Y b Y coordinate for the b ,h boring sample (Appendix A-84) 

d^ Available moisture depth of the soil mass (Appendix A-90) 

d hbar Size of horizontal barrier (Appendix A-86) 

^source Distance effect is away from the edge of barrier/slab (Appendix A-86) 

d^ Depth of vertical barrier (Appendix A-86) 

d zone Depth of the root or flower bed zone (Appendix A-86) 

h c c* suction measurement in a sample (Appendix A-87) 

h dry Dry suction limit for soils in the field, bare soil at the surface usually 

has a value of 6.0pF and soil with vegetation at the surface usually 
has a value of 4.5pF (It may be requested in various algorithms. 
Appendix A-90 is the first time h dry is requested as input). 
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h^i Wet limit suction for soils in the field, usually 2.5pF (It may be requested 

in various algorithms. Appendix A-90 is the first time h w , is used), 
n Local annual weather cycle (Appendix A-83). 

t const The month the slab was constructed (Appendix A- 1 07) 

t dry Driest month of the year (Appendix A- 1 07) 

t^ The month the soil volume change is sought (Appendix A- 1 07) 

p s Percent passing 2 micron for layer s (Appendix A-84) 

v s Percent passing number 200 sieve for layer s (Appendix A-84) 

More detailed explanations of how these variables are used is described in the following 
sections of this document. Appendix A gives a complete process flow detailing the use 
of the variables listed above and the variables listed in the list of symbols. 

Input, as requested by the process flow in Appendix A and as graphically 
represented in the software design layout sheets in Appendix B, is broken down into six 
major categories: 

1. Project Information input. This includes project name, date the project started, 
name of the engineer assigned to the project, and the site location. This information 
is basic information an engineering firm may want to establish in a database so that 
information can be readily accessed. The only information critical to calculating 
differential soil movement is the site location. The site location can be used as the 
means to locate the TMI for a site without requiring the user to manually research 
this bit of information. Refer to Appendix A-82 for a process flow and Appendix B- 
117 for a program design layout. The TMI for a particular site is found in Fig.26. 

2. Foundation Layout. This is basic surveying information to establish a three 
dimensional grid for logging information related to appropriate physical locations. 
Refer to Appendix A-83 for process flow and Appendix B-l 1 8 for program design 
layout. 

3. Soil Layer Property Input. These values include boring sample data and soil layer 
property data. The data entered here will provide enough information to generate a 
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three dimensional picture of the soil profile. Refer to Appendix A-84 for a process 
flow and Appendix B-l 19, 120, 126 for program design layouts. The soil profile 
used in all example calculations is shown in Table 1. 




Fig. 26. (Wray 1978) Thornthwaite Index for Texas (20 year average 1955-1974) 



Table 1. Soil Profile Considered in all Sample Calculations 



Soil 


L s 


LL S 


PL S 


V s 


Ps 


Yd 


Layer No. 


ft 


% 


% 


% 


% 


lbs/ft 3 


s = 1 


5 


56.0% 


22.0% 


65.0% 


25.0% 


120 


CM 

II 

(/) 


9 


67.0% 


24.0% 


76.0% 


46.0% 


120 


s — 3 


13 


44.0% 


18.0% 


47.0% 


18.0% 


120 


s = 4 


30 


53.0% 


20.0% 


64.0% 


31.0% 


120 
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4. Moisture Effect Input. This gives the user the option to input the type and location 
of the moisture effect cases considered in vertical differential soil movement 
calculations. Refer to the plan view, Appendix C-l 32, for a view of the foundation 
used as the example throughout this document. The input for location and moisture 
effect case type can be found listed on the program design layout Appendix B-121 . 
For example, case “1-5” describes the sidewalk with grass covered lawn on the 
north and west portions of the building shown in Appendix C-l 32. The number “1”, 
in “1-5”, describes this case as the first sequential case of this type. The number 
“5”, in “1-5”, describes this case to be grass at the surface with a horizontal barrier. 
Additional information is shown in Appendix B-121 and is described in the process 
flow sheets Appendix A-85, 86. 

5. Measured Suction Profile Input. Geotechnical engineers may identify some 
unique moisture effect cases that cannot be modeled in the procedures to address the 
typical moisture effect cases contained herein. Those special cases are addressed 
using initial and final measured suction profiles. The input of these measured 
suction profiles is described in the process flow sheet Appendix A-87, and shown in 
the program design layout sheet, Appendix B-122, 123. 

6. Time Input for Post Construction Cases. Algorithms for post construction 
solutions are time dependent. Calculating the post construction differential soil 
movements can only be accomplished if the driest (or wettest) time of year is 
known, t dry , the month of year the foundation was built is known, t const , and the month 
the soil volume change is sought is known, t vc . Input of these quantities is described 
in the process flow sheet Appendix A- 107. 

Calculated Soil Properties and Soil Characteristics 

Once the liquid limit, the plastic limit, the percent passing the number 200 sieve, 
the percent smaller than two microns, and the dry unit weight is entered for the soil 
profile considered, the remainder of the required soil properties and characteristics are 
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solved. The process flow sheet, Appendix A-88, describes the algorithm by which these 
remaining soil characteristics are solved. Definitions for all symbols are found in the 
list of symbols. 

The plasticity index, PI(Zj), is the difference between the liquid limit, LL(Zj), and 
the plastic limit, PL(Zj). The percent fine clay, Clay(Zj), is calculated by dividing the 
percent smaller than two microns, p(Zj), by the percent passing the number 200 sieve, 
v(Zj). The activity ratio, Ac(Zj), is calculated by dividing the PI(Zj) by Clay(z,). The 
cation exchange capacity, CEC(Zj), is the larger of the values PL(Zj) 1 17 or LL(Zj) 0912 . 

The cation exchange activity, CEAc(Zj), is calculated by dividing CEC(Zj) by Clay(z,). 
Using Fig. 15, the matrix suction compression index, y hloo , is the value found at the 
intersection of CEAc(Zj) and Ac(Zj). To obtain the y h (Zj) simply multiply Clay(zi) by 
y hl00 . The slope of the suction versus the gravimetric water content line, S(Zj), is 
approximated by (16). Mitchell’s diffusion coefficient, a(Zj), is approximated by (17). 

Mitchell’s unsaturated permeability, K 0 (Zj) |h 0 | , is calculated by (18). Using Fig. 16 
for center lift conditions and Fig. 17 for edge lift conditions, the edge moisture variation 
is found by finding the curve that corresponds with the soil's diffusion coefficient, a(Zj), 
then finding the intersection made by drawing a vertical line corresponding to the local 
TMI (for locations in Texas the local TMI can be found in Fig. 26), and finally drawing 
a horizontal line from that intersection across to the y axis identifies the appropriate 
edge moisture variation distance. For each case of center lift and edge lift, the largest 
edge moisture variation value, considering all incremental soil layers, is the soil 
profile’s edge moisture variation distance, e mc and e me . 

Procedures to Calculate the Equilibrium Suction 

The next step in the process of calculating the vertical differential soil movement 
is to establish a baseline equilibrium suction, h,,, (units of pF). All suction profiles are 
calculated relative to this baseline. The equilibrium suction, h,,,, is calculated for the soil 
layer located at the depth to constant suction. The depth to constant suction, z is 
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dependent upon the h m Therefore, if the depth to constant suction, z,,,, does not fall 
within the soil layer in which the h m was calculated, then the h in must be re-calculated 
for the soil layer located at the depth to constant suction. This is an iterative process. 
The process flow sheet which describes the procedure for calculating h, n is found on 
Appendix A-90. The procedures for calculating z*, are found in the next section. The 
equilibrium suction value is dependent upon the following inputs: TMI, d ain , h wl , h diy , 0 S , 
0 r , A, and B. For the state of Texas the TMI is obtained from Fig. 26. Values for 0 S , 0 r , 
A, and B for a particular classification of soil are listed in Table 2. As reported by Gay 
(1994), a conservative estimate for the d am can be taken as 30 cm for all expansive soils 
in which the d ain is not known. All variables have default values. Refer to Appendix A- 
90. Once all input is entered y, d„ T,, 0 W „ 0 dry are calculated using (25), (26), (27), (38), 
and (39) respectively. These values are used in equations (33) and (37) to solve for h,,,. 

A value of is picked and an iterative process for solving for 0 m is performed until 
equal values of 0 m satisfies both equations (33) and (37). 

Depth to Constant Suction 

The depth to constant suction is the depth in a soil profile to which there is no 
longer a significant seasonal suction change. The definition of no significant seasonal 
suction change is 0.2 pF from Lytton’s 1994 paper. Depth to constant suction is 
dependent upon the diffusivity of the soil, the number of local annual weather cycles, 
and the suction amplitude at the surface. As discussed in the “Soil Volume Strain” 
section, the soil’s ability to change volume is dependent upon the seasonal change of 
moisture and the magnitude of the mean principal stress present to suppress the volume 
change. Using Mitchell’s equation (42) there will always be a reduction in the suction 
amplitude with an increase of depth. Additionally, an increase in depth introduces 
added restraint to soil volume strain due to the effect of overburden. Considering these 
two aspects of the soil volume change equation it becomes apparent that for depths 
greater than 30 feet there is usually no longer any significant changes in suction. The 
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depth to constant suction, z will be used as the depth beyond which there is no 
significant change in suction. 

Table 2. (Lytton et al. 1989) Gardner’s Coefficients 



No 


Unified 

Soil 

Class. 


Void 

ratio 

e 


AWL 


XWL 


Porosity 

n = e s 


0 r =O.10 s 


A 


B 


GM-GC1 


GM-GC 


0.42 


0.004 


0.637 


0.296 


0.0296 


1083.777 


0.637 


GM1 


GM 


0.259 


0.152 


0.269 


0.206 


0.0206 


12.2224 


0.269 


GM2 


GM 


0.382 


0.04 


0.648 


0.276 


0.0276 


111.1578 


0.648 


GM3 


GM 


0.451 


0.066 


0.251 


0.311 


0.0311 


27.00574 


0.251 


GM4 


GM 


0.608 


0.043 


0.478 


0.378 


0.0378 


69.90875 


0.478 


GP1 


GP 


0.255 


0.065 


0.55 


0.203 


0.0203 


54.58668 


0.55 


GW1 


GW 


0.443 


0.039 


0.302 


0.307 


0.0307 


51.39672 


0.302 


GW2 


GW 


0.506 


0.596 


0.318 


0.416 


0.0416 


3.489424 


0.318 


GW3 


GW 


0.252 


0.309 


0.319 


0.355 


0.0355 


6.745925 


0.319 


SM-SC1 


SM-SC 


0.656 


0.013 


0.77 


0.396 


0.0396 


452.9567 


0.77 


SMI 


SM 


0.396 


0.016 


0.562 


0.284 


0.0284 


227.9712 


0.562 


SM2 


SM 


0.619 


0.001 


1.023 


0.382 


0.0382 


10543.87 


1.023 


SM3 


SM 


1.192 


0.11 


0.339 


0.544 


0.0544 


19.843 


0.339 


SM4 


SM 


0.655 


0.039 


0.468 


0.396 


0.0396 


75.32435 


0.468 


SM5 


SM 


0.623 


0.015 


0.835 


0.384 


0.0384 


455.9411 


0.835 


SM6 


SM 


1.024 


0.011 


0.671 


0.506 


0.0506 


426.194 


0.671 


SM7 


SM 


0.72 


0.023 


0.549 


0.419 


0.0419 


153.9119 


0.549 


SM8 


SM 


0.831 


0.208 


0.436 


0.454 


0.0454 


13.12009 


0.436 


SM9 


SM 


0.721 


0.01 


0.835 


0.419 


0.0419 


683.9116 


0.835 


SM10 


SM 


0.818 


0.018 


0.806 


0.45 


0.045 


355.4082 


0.806 


SM1 1 


SM 


1.132 


0.029 


0.745 


0.531 


0.0531 


191.6911 


0.745 


SM12 


SM 


2.16 


0.042 


0.501 


0.684 


0.0684 


75.46589 


0.501 


SP-SM1 


SP-SM 


0.574 


0.048 


0.769 


0.365 


0.0365 


122.3936 


0.769 


SP-SM2 


SP-SM 


0.583 


0.095 


0.613 


0.368 


0.0368 


43.17938 


0.613 


SP1 


SP 


0.712 


0.0554 


0.79 


0.416 


0.0416 


111.2987 


0.79 


SP2 


SP 


0.736 


0.024 


0.951 


0.424 


0.0424 


372.2106 


0.951 


SP3 


SP 


0.818 


0.042 


0.9 


0.45 


0.045 


189.1258 


0.9 


SP4 


SP 


0.593 


0.076 


0.665 


0.372 


0.0372 


60.83961 


0.665 


SP5 


SP 


0.583 


0.053 


0.809 


0.368 


0.0368 


121.5414 


0.809 


SW-SP1 


SW-SP 


0.773 


0.012 


1.082 


0.296 


0.0296 


1006.512 


1.082 


SW1 


SW 


0.551 


0.162 


0.5852 


0.355 


0.0355 


23.75117 


0.5852 


SW3 


SW 


0.385 


0.111 


0.616 


0.278 


0.0278 


37.21149 


0.616 


CL1 


CL 


1.183 


0.065 


0.417 


0.542 


0.0542 


40.1871 


0.417 


CL2 


CL 


0.621 


0 


0.976 


0.383 


0.0383 


0 


0.976 


CL3 


CL 


0.785 


0.004 


0.61 


0.44 


0.044 


1018.451 


0.61 


CL4 


CL 


0.61 


0 


0.957 


0.379 


0.0379 


0 


0.957 


CL5 


CL 


0.825 


0.018 


0.523 


0.452 


0.0452 


185.2369 


0.523 


CL6 


CL 


0.964 


0 


1.91 


0.491 


0.0491 


0 


1.91 


CL7 


CL 


0.808 


0 


1.593 


0.447 


0.0447 


0 


1.593 


CL8 


CL 


0.972 


0.028 


0.344 


0.493 


0.0493 


78.85731 


0.344 
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Table 2. Continued 



No 


Unified 

Soil 

Class. 


Void 

ratio 

e 


AWL 


XWL 


Porosity 

n = 0 S 


0=0.10* 


A 


B 


CL9 


CL 


0.562 


6.9E-06 


1.551 


0.36 


0.036 


5154077 


1.551 


CLIO 


CL 


0.517 


0.024 


0.359 


0.341 


0.0341 


95.23328 


0.359 


CL11 


CL 


0.415 


0.064 


0.247 


0.293 


0.0293 


27.59434 


0.247 


ML8 


ML 


0.858 


0.003 


0.973 


0.462 


0.0462 


3132.411 


0.973 


ML9 


ML 


0.606 


0.022 


0.707 


0.377 


0.0377 


231.514 


0.707 


ML10 


ML 


2.163 


0.038 


0.61 


0.684 


0.0684 


107.2053 


0.61 


ML1 1 


ML 


0.603 


0.031 


0.722 


0.376 


0.0376 


170.0741 


0.722 


ML12 


ML 


0.821 


0 


2.059 


0.451 


0.0451 


0 


2.059 


CL12 


CL 


0.972 


0.004 


0.474 


0.493 


0.0493 


744.6291 


0.474 


ML-OL1 


ML-OL 


1.721 


0.064 


0.535 


0.632 


0.0632 


53.55747 


0.535 


ML-OL2 


ML-OL 


0.885 


0 


1.362 


0.469 


0.0469 


0 


1.362 


ML-CL1 


ML-CL 


0.692 


0.066 


0.365 


0.409 


0.0409 


35.11204 


0.365 


ML-CL2 


ML-CL 


0.642 


0.013 


0.634 


0.391 


0.0391 


331.1743 


0.634 


ML-CL3 


ML-CL 


0.645 


0.001 


1.052 


0.392 


0.0392 


11271.97 


1.052 


ML1 


ML 


1.015 


0.003 


1.054 


0.504 


0.0504 


3774.668 


1.054 


ML2 


ML 


1.067 


0 


1.257 


0.516 


0.0516 


0 


1.257 


ML3 


ML 


1.119 


0.032 


0.681 


0.528 


0.0528 


149.9167 


0.681 


ML4 


ML 


1.558 


0.065 


0.411 


0.609 


0.0609 


39.63571 


0.411 


ML5 


ML 


0.88 


0 


1.019 


0.468 


0.0468 


0 


1.019 


ML6 


ML 


0.825 


0 


1.474 


0.452 


0.0452 


0 


1.474 


ML7 


ML 


0.783 


0.012 


0.575 


0.439 


0.0439 


313.1978 


0.575 



The depth to constant suction varies dependent upon the soil strata, the local 
weather conditions, and the type of soil surface conditions. The depth to constant 
suction, z^, will be calculated for the following cases of soil surface conditions, 
Appendix A-91 : 

1 . Bare soil at the surface (Appendix A-92). 

2. Grass at the surface with shallow roots (Appendix A-92). 

3. Flower bed at the surface with a known flower bed zone depth, d zone (Appendix A- 
93). 

4. A tree at the surface with a known depth of root zone, d zone (Appendix A-93). 

5. And measured suction profiles (Appendix A-94). 

The simplest case would be that of bare soil or grass at the surface with only one 
soil layer the entire depth of moisture active zone. For this case equation (46) is all that 
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is needed to solve for 2 ^. Simply replace the numerator in the natural log function with 

0. 2. pF, substitute entered or previously calculated values for all other variables and 
solve for z^ The suction amplitude, U 0 , is the difference between maximum and the 
minimum surface soil suction values. For bare soils the suction amplitude is U 0 = h d[y - 
h^ where h dry = 6.0 pF. For soils with grass at the surface the suction amplitude is U 0 = 
h^ - h^, where h dry = 4.5 pF. Refer to the “Typical Suction Levels in the Field” for 
typical values of suction. 

In most soil profiles there is usually at least two soil layers to consider. In cases 
considering bare soil or grass at the surface, with multiple soil layers, the algorithms in 
Appendix A-92 are used to solve for z, n . The procedures are simple. 

1 . First solve for the ratio of the suction difference between the top increment and 
bottom increment of a soil layer, P s , using equation (45). 

2. Multiply P s by the suction amplitude at the top of the layer which gives the suction 
amplitude at the bottom of the layer, U Ls . 

3. Test U Ls : 

a. If U Ls is less than 0.2 pF, then the depth to constant suction is in this 
layer, in which case z ^ is solved using (46). The initial U 0 is the same as 
the one described in the previous paragraph. All subsequent U 0 ’s are 
solved using the procedure in paragraph 3.b. described below. 

b. If is grater than 0.2 pF, then set U 0 equal to and repeat these three 
steps. 

The depth to constant suction is the z,,, calculated for the last layer added to the value of 
the depth to the top of this layer. 

The procedures for calculating the z, n for flower beds and trees is basically the 
same with two distinct differences: (1) The suction amplitude is the difference between 
the limit suction value and the equilibrium suction value. (2) The suction amplitude at 
the surface is constant to the depth of moisture effect zone, d zone . For a flower bed at 
the surface the suction amplitude is the difference between the equilibrium suction value 
and the wet limit suction value, U 0 = !!„, - h w] . For a tree at the surface the suction 
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amplitude is the difference between the dry limit suction value for vegetation and 
equilibrium suction value, U 0 = h dry - h,,,. Refer to Appendix A-93 for a detailed process 
flow of these procedures. Refer to “Typical Suction Levels in the Field” and 
“Procedures to Calculate the Equilibrium Suction” for values of (h^, and h dry ) and h^ 
respectively. 

When making calculations of the depth to constant suction, z^, for the case of 
measured suction profiles, the depth to constant suction can be taken as the value of z ^ 
calculated for one of the four previously explained cases if the surface conditions are the 
same. However, a check must be made to ensure that the z^ in a particular measured 
suction profile case does not extend beyond those values of z, n in the previously 
described cases. If the z, n extends beyond those values calculated using the procedures 
described in the previous four cases then the z ^ must be calculated using the procedures 
described in Appendix A-94. The procedures are the same as those described above for 
bare soil and grass at the surface with one major difference. The difference is that the 
process is initiated at the bottom of the measured suction profile. This z ^ value 
calculated for a measured suction profile is then compared to the z ^ for the surface 
conditions most similar to the measured profile surface conditions. The depth to 
constant suction is taken as the larger of these values. 
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POST EQUILIBRIUM AND POST CONSTRUCTION SOLUTIONS 

At this point in the procedures all inputs and calculated values required to generate 
suction profiles for soil volume strain calculations have been explained. The procedures 
to generate soil suction profiles and the corresponding soil movement for moisture effect 
cases and design effects are explained using the following two procedural approaches: (1) 
the post equilibrium approach, and (2) the post construction approach. The post 
equilibrium approach refers to procedures to solve for differential movement of expansive 
soils after the soil under the slab has reached an equilibrium moisture content. When a 
slab is placed on an expansive soil, there is a time period in which the moisture content 
under the slab changes from the suction profile as it was when the slab was placed to 
approach the equilibrium suction value. The algorithms to solve for differential soil 
movement during this time period are handled by the post construction approach. 

Post Equilibrium Approach 

Once the soil under the slab reaches the equilibrium suction value, the differential 
soil movement is calculated using the differential in soil suction profiles between the 
equilibrium suction profile under the slab and the suction profile at the edge of the slab. 
The suction profile at the edge of the slab is dependent upon the seasonal moisture change. 
The restraint to soil volume increase at any depth is dictated by the amount of overburden 
and surcharge stress. The post equilibrium algorithms to solve for the limit cases are: (1) 
differential soil movement causing a maximum center lift distortion mode, and (2) 
differential soil movement causing a maximum edge lift distortion mode. 

Equilibrium Suction Profiles 

The first calculations that are made in the post equilibrium approach are the base 
line equilibrium suction profiles. The equilibrium suction profile is generated by reducing 
the equilibrium suction value calculated for a specific moisture effect case by one 
centimeter in suction for every centimeter of descent in soil depth. Care should be taken to 
convert the calculated equilibrium suction value from units of pF to units of centimeters. 
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Once the equilibrium suction profile is generated these values of suction in units of 
centimeters should be converted back to units of pF to make further calculations. Refer to 
Appendix A-96 for algorithms to calculate equilibrium suction profiles. 

Suction Profiles Causing The Center Lift Distortion Mode 

Center lift distortion modes are caused by a differential in the equilibrium suction 
profile under the slab and the dry limit suction profile at the edge of the slab. Once the 
equilibrium suction profiles are calculated for each case considered, the dry limit suction 
profiles for the slab edge moisture effect conditions are calculated. The dry limit suction 
profiles for bare soil and grass at the slab edge are calculated using Appendix A-97, using 
variables unique to the specific moisture effect case and location. The dry limit suction 
profiles for a flower bed and a tree at the edge of slab are calculated using the Appendix A- 
98. The suction profiles are calculated by the application of (43). The differences in these 
processes, dependent upon moisture effect case and location are: the surface suction values, 
the soil properties, the baseline equilibrium suction profiles, and depth to constant suction 
values. For all moisture effect cases, the suction amplitude change and the depth are 
initialized at every soil layer change. The cases of a tree and a flower bed at the surface 
have one other major difference, these cases include the depth of moisture effect zone. The 
suction profile within the tree root or flower bed zone is constant and equal to the limit 
suction value for the full depth of the moisture effect zone, d zone . Refer to Fig. 27 for an 
example of the dry limit post equilibrium suction profiles for a flower bed at the surface 
with a 4 foot moisture effect zone depth. Refer to Fig. 28 for an example of the dry limit 
post equilibrium suction profiles for a tree at the surface with 20 foot deep roots. 

Suction Profiles Causing The Edge Lift Distortion Mode 

Edge lift distortion modes are caused by a differential in the equilibrium suction 
profile under the slab and the wet limit suction profile at the edge of the slab. Once the 
equilibrium suction profiles are calculated for each case considered, the suction profiles for 
the slab edge moisture effect conditions are calculated. The wet limit suction profiles for 
bare soil, grass, and trees at the slab’s edge are calculated using Appendix A-99. The wet 
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limit suction profiles for a flower bed at the edge of the slab is calculated using the 
algorithm described in Appendix A- 100. The differences in these processes which depend 



Suction (pF) 

IO IO O « 4* ui 




Fig. 27. Post Equilibrium Suction Profiles for Flower Bed at the Surface (Dry Limit) 
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Fig. 28. Post Equilibrium Suction Profiles for a Tree at The Surface (Dry Limit) 
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upon the moisture effect case and location are: the soil properties, the baseline equilibrium 
suction profiles, and depth to constant suction values. Note that the suction profile is 
constant and equal to the limit suction value for the full depth of the flower bed zone, d zone . 
The suction amplitude change and the depth are initialized at every soil layer change. 
Refer to Fig. 29 for an example of the wet limit post equilibrium suction profiles for a 
flower bed at the surface. 
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Fig. 29. Post Equilibrium Suction Profiles for a Flower Bed at The Surface (Wet Limit) 



Vertical Differential Soil Movement 

Expansive soil volume change is predicted by applying the theory explained in the 
“Soil Volume Strain” section of the previous chapter. The moisture effect cases considered 
consist of: 

1 . Bare soil at the surface, 

2. Grass at the surface with shallow roots, 
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3. A flower bed at the surface with a known depth of flower bed zone, d zone , and 

4. A tree at the surface with a known depth of root zone, d zone 

The procedures also consider the application of design effects to provide the best possible 
conditions so foundation design performance is at a pinnacle. The design effects 
considered include: 

1. Vertical barriers, 

2. Horizontal barriers, and 

3. The effect of locating the moisture effect case a distinct distance from the slab edge. 
Explanation of the procedures to predict vertical differential soil movement is broken down 
into three categories: (1) all moisture effect cases with or without a vertical barrier, (2) 
moisture effect cases with or without a horizontal barrier, and (3) moisture effect cases 
located a distinct distance from the slab edge. 

Post equilibrium vertical differential soil movement can be determined for all 
moisture effect cases with or without a vertical barrier by simply applying the algorithms 
described by Appendix A- 103. This process can be broken down into three basic steps: 

(1) All pertinent variables and information collected or calculated in previous procedural 
steps are retrieved for use in this process. (2) Tests are performed on variables to ensure 
the appropriate calculations and quantities are used, based on Lytton’s volume strain 
theory. (3) Incremental strain and corresponding vertical soil movement is summed up for 
every soil increment within the moisture active zone of the soil profile. The total vertical 
differential soil movement quantity is defined by the variable Az. Refer to Fig. 30 for an 
example of bare soil dry limit suction profiles (center lift distortion mode) for a slab edge 
that has a four foot deep vertical barrier. 

Post equilibrium vertical differential soil movement can be computed for all 
moisture effect cases, except the tree case, with or without a horizontal barrier by simply 
applying the algorithms described by Appendix A- 104. The case of the tree is unique in 
that the tree roots will be active regardless of whether a horizontal barrier exists or not. If 
the tree is killed because of the placement of the horizontal barrier then the point becomes 
mute. The process in Appendix A-104 is broken down into four basic steps: (1) All 
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Fig. 30 Bare Soil Dry Limit Suction Profiles for a Slab with a 4ft Vertical Barrier 

pertinent variables and information collected or calculated in previous procedural steps 
are retrieved for use in this process. (2) The horizontal velocity for each incremental soil 
layer is calculated based on the edge moisture variation distance and the size of horizontal 
barrier. This step is based on equation (52). (3) The suction at the slab edge is calculated 
using equation (53). (4) This slab edge suction value, h edge (z,), is used in the place of the 
final suction value in Appendix A- 103 to make vertical differential soil movement 
calculations, Az. Refer to Fig. 3 1 for an example of wet limit suction profiles (edge lift 
distortion mode) for a four foot deep flower bed at the edge of a four foot horizontal barrier 
constructed adjacent to the edge of the slab. 

If a flower bed or tree is located a distinct distance from the slab edge, the vertical 
differential soil movement is calculated by applying Appendix A- 104. These are the only 
two cases that can be located a distance from the slab and still have an effect on vertical 
differential soil movement under the slab. The differential movement caused by the tree or 
flower bed located a distance from the slab edge should be compared to the differential 
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Fig. 31. Suction Profiles for Flower Bed at the Slab Edge with a 4 ft Horizontal Barrier 

movement caused by the surface condition that exists between the slab edge and the 
tree/flower bed. The maximum differential soil movement produced by either the 
tree/flower bed case or the conditions between the slab edge and the tree/flower bed case 
controls. This process is broken down into four basic steps: (1) All pertinent variables and 
information collected or calculated in previous procedural steps are retrieved for use in this 
process. (2) The horizontal velocity for each incremental soil layer is calculated based on 
the edge moisture variation distance and the distance the moisture effect is away from the 
slab edge, d source This step is based on equation (52). (3) The suction at the slab edge is 
calculated using equation (53). (4) This slab edge suction value, h edge (Zj), is used in the 
place of the final suction value in Appendix A- 103 to make vertical differential soil 
movement calculations, Az. Refer to Fig. 32 for a dry limit suction profiles (center lift 
distortion mode) for tree roots that are 20 ft deep and located four feet away from the edge 
of a slab with a four foot vertical barrier. 
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Fig. 32. Suction Profiles for a Tree 4ft From the Slab Edge with a 4ft Vertical Barrier 
Post Construction Approach 

As described in the “Post Construction Theory” section earlier, the worst case of 
vertical soil movement in expansive soils may occur in the post construction time period of 
the life of a structure. Lytton and Woodbum’s 1973 paper provides a perfect illustration of 
this worst case scenario occurring in the post construction time period during the life of a 
school building. The procedures to apply the post construction theory is basically a two 
step process. First, all suction profiles considered are calculated. Second, these suction 
profiles are used to calculate vertical soil movement at the center of the slab, Az slab , and 
combined with the vertical soil movement at the slab edge, Az edge , to arrive at a total 
vertical differential soil movement quantity, Az, otal . 

Post Construction Suction Profiles 

The suction profiles must be calculated for two time periods: soil suction profiles 
at the time of construction, t cons „ and soil suction profiles at the time the vertical differential 
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soil movement is sought, t vc . Soil suction profiles at the edge of the slab are generated for 
the time at which the vertical differential soil movement is sought. Soil suction profiles 
beneath the slab are generated at the time of construction and at the time the vertical 
differential soil movement is sought. 

Appendix A- 107 is the process flow sheet that: describes the time variables, lists a 
step by step process for entering the required times, and lists the steps to calculate the time 
periods needed to make soil suction profile calculations. Appendix B-127 is a picture of 
the program design layout sheet for the time input and display of suction profiles both 
under the slab and at the edge of the slab for time that transpires from the driest month to 
the month the volume change is sought. In Appendix B-127, the suction profile under the 
slab at the time of construction is the thin solid line. The dashed line is the suction profile 
under the slab six months after construction. Notice that this dashed line is tending toward 
equilibrium. The equilibrium profiles are the vertical thin lines with tick marks. Finally, 
the thick solid line is the suction at the edge of the slab with a vertical barrier. To generate 
suction profiles illustrated in Appendix B-127 the following time inputs are required: (1) 
the driest month of the year, t^, (2) the month the slab is constructed, t const , and (3) the 
month the soil volume change is sought, t vc . These entered values of time are required for 
two primary purposes: (1) to establish a baseline time to enable application of the surface 
annual weather cycle model developed by Mitchell, (40), and (2) to establish time periods 
used in the procedures to calculate suction profiles for each case at times required to make 
calculations for the post construction algorithms. Appendix A- 107 is broken down into two 
basic steps. First, the times required for input are entered. Second, three time periods are 
calculated using the driest month of the year as the baseline. The time period used to 
calculate the suction profile under the slab at the time of construction is represented by the 
variable t cons(/dry , as calculated in Appendix A- 107. The time period between the time of 
construction and the time at which the soil volume change is sought is represented by the 
variable t^,^,, as calculated in Appendix A- 107. Finally, the variable t^^ describes the 
time transpired from the driest month of the year to the month the soil volume change is 
sought. 



67 



All soil suction profiles with respect to surface conditions, time and depth are 
calculated using the algorithms in Appendix A- 109 and Appendix A-l 10. The process by 
which post construction soil suction profiles are generated consist of the following basic 
steps: 

1 . All values required to make calculations are retrieved for each particular moisture 
effect case and time the suction profile is sought, Appendix A- 109. 

2. Suction values, for each particular moisture effect case, are calculated for 5 centimeter 
thick incremental soil layers. 

3. For every incremental soil layer cycle, variables are re-initialized and pertinent 
equations are activated based on: the location within the soil profile, the moisture effect 
case type, and the type and size of design effect implemented (If a horizontal barrier is 
implemented or if a moisture effect case is located a distance from the edge, Appendix 
A- 104 will be used to augment the procedures described in Appendix A-l 10.) 

Once all soil suction profiles are generated, the soil suction profile at the slab center 
is damped from time t const to time t vc by applying the algorithms described by Appendix A- 
111. This process takes the soil suction profile under the slab at the time of construction, 
h slab (Zj,t const/diy ), and dampens it based on equations (59) and (61). These equations model 
the soil suction profile under the center of the slab at time t const so that the soil suction in 
each incremental soil layer approaches the equilibrium suction profile. The rate of 
approach depends on unsaturated soil profile properties, surface weather conditions, and 
time. The process in Appendix A-l 1 1 is broken down into two basic steps. First, all 
required variables are retrieved. Second, the dampened suction value for every incremental 
soil layer is calculated for the time transpired from the time of construction to the time the 
soil volume change is sought. 

Vertical Differential Soil Movement 

Once all soil suction profiles are calculated, the expansive volume strain theory 
developed by Lytton is applied. All moisture effect cases are applied to the same vertical 
differential soil movement algorithms which are broken down into a four step process. 
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First, all variables and quantities required for calculations are retrieved, as described in 
Appendix A-l 13, depending on which moisture effect case and design effect scenario is 
chosen. Second, vertical differential soil movement is calculated for soil at the edge of the 
foundation by way of steps described in Appendix A-l 14. Third, vertical differential soil 
movement is calculated for soil under the center of the foundation by way of steps 
described in Appendix A-l 15. Finally, the two quantities of vertical differential soil 
movement are summed together to arrive at the total differential soil movement for the soil 
under the slab at a specified time after construction, as described in Appendix A-l 15. 

Appendix B-128 is a program design layout sheet that illustrates the results of the 
post construction procedures when applied to calculate the worst case center lift mode for 
case “1-6”. The moisture effect case considered in Appendix B-128 is grass at the surface 
with a 4 ft vertical barrier, case “1-6”. The vertical soil movement calculated at the edge 
moisture variation distance under the slab is Az slab = -0.13 inches (shrink). The vertical soil 
movement calculated at the edge of the slab is Az edge = -0.61 inches. The total vertical 
differential soil movement for this moisture effect case using Table 2 as the soil profile is 
Az, otal = 0.48 inches (shrink). This quantity, Az, ota | , is the difference between the values of 
Az slab and Az edge . 
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MEASURED SUCTION PROFILES 

When cases exist that can not be modeled by one of the moisture effect cases 
discussed in the previous procedures, vertical differential soil movement can be calculated 
by the measured suction profile method. This is a method to solve for the differential soil 
movement by measuring the initial and final suction profiles and applying those profiles to 
the volume strain theory developed by Lytton. The procedures are broken down into three 
broad steps: (1) The measured suction values for a sample are entered and stored for later 
use, Appendix A-87. (2) The entered measured suction values for a sample are converted 
into a measured suction profile, Appendix A- 101. (3) The measured suction profiles are 
used to solve for vertical differential soil movement for a specific column of soil, 

Appendix A- 105 and Appendix A- 106. 

The first step to applying this methodology is to input measured suction values for 
a sample of soil. There are a number of methods to measure the matrix suction in a soil 
sample. One of the more accurate means for measuring total soil suction is the 
thermocouple psychrometer. The thermocouple psychrometer sensor and data logging 
equipment is not inexpensive and requires some expertise and maintenance to operate. For 
these reasons the filter paper is becoming a more widely accepted method in practice. The 
filter paper is inexpensive, requires practically no maintenance, and requires very little 
expertise to obtain results nearly as accurate as the thermocouple psychrometer (McKeen 
1981). 

The process of inputting measured suction values occurs as part of the process to 
input moisture effect cases. Appendix B-121 is a program design layout sheet that 
illustrates the choices for the moisture effect case input. The program user would select 
“12. Measured Suction Profile” and the next screen that appears is illustrated by Appendix 
B-122. Appendix B-122 and 123 are a program design layout sheets that illustrate 
windows environment screens for inputting the measured suction profiles. Appendix A-87 
is the process flow describing measured suction profile input. It is important to note that a 
surface suction value must be entered or a default value equal to the equilibrium suction 
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will be assumed. 

Appendix A- 101 is the process flow describing the steps to generate a measured 
suction profile based on the entered suction values. The procedures are broken down into 
two basic steps: (1) The slope of a line connecting adjacent suction values within a 
sample is calculated. (2) The measured suction profile is generated, using this previously 
calculated slope, by connecting a line to adjacent suction values in the sample. Table 3 is 
an example of values used to generate a measured suction profile. Fig. 33 is a graph of the 
initial and final measured suction profiles generated by the values listed in Table 3 



Table 3 Measured Suction Profile Values 



Input Suction Values 




Initial 






Final 






(ft) 

Depth 


(cm) 

Depth 


(PF) 

Suction 


(cm) 

Suction 


(ft) 

Depth 


(cm) 

Depth 


(PF) 

Suction 


(cm) 

Suction 


0 


0 


4.5 


31622.78 


0 


0 


2.5 


316.2278 


3 


91.44 


4.3 


19952.62 


3 


91.44 


3 


1000 


12 


365.76 


4 


10000 


12 


365.76 


3.3 


1995.262 


16 


487.68 


3.7 


5011.872 


16 


487.68 


3.5 3162.278 


20 


609.6 


3.57 


3715.352 


20 


609.6 


3.57 


3715.352| 



Slope 

Initial Final 



(ft) 

Depth 



°l 


-0.00219 


0.005468 


3 | 


-0.00109 


0.001094 


12 | 


-0.00246 


0.00164 


16 | 


-0.00107 


0.000574 


20 



Given Values 


z m = 


20 ft 


“ 


3.57 pF 


Gi = 


40 cm 


K 0 (dry) ” 


0.33 


K 0 (wet) ” 


0.67 


f<dry) = 


0.5 


f (wet) ™ 


0.8 



Applying the measured suction profiles to predict vertical differential soil 
movement is probably the simplest of all procedures contained herein. The procedures are 
broken down into two basic steps: (1) select an initial and final suction profile to 
accurately model the moisture change scenario being considered. Appendix A-105, and (2) 






apply the algorithms to calculate vertical differential soil movement based on these 
measured suction profiles, Appendix A- 106. 
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Fig. 33. Initial and Final Measured Suction Profiles Described in Table 4 



The challenge is in deciding what suction profiles are used as the initial and final 
suction profile to accurately model the soil conditions which produce reasonable limit 
conditions used for design. The answer should come from knowledge about history of the 
site, to include: surface drainage conditions, soil characteristics in the soil profile, surface 
weather conditions, and type of structure placed. For example, consider the construction 
site described by the following: the site is located in College Station, the site is surface 
drained, the soil profile is described by Table 1, and the measured suction profile is located 
where a large oak tree was removed as a result of pre-construction clearing. The measured 
suction profile at this location is a dry profile. This profile will tend toward equilibrium 
and may approach the wet limit state subsequent to a wet season. A conservative approach 
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for calculating vertical differential soil movement is to set the initial suction profile equal 
to the measured suction profile and set the final suction profile equal to the wet limit state 
for the surface conditions that exist after the construction project is complete. This 
scenario produces heave values of vertical differential soil movement. Table 3 and Fig. 33 
are values used to obtain results for vertical differential soil movement in the example 
described above. The vertical differential soil movement in this case is 5.48 inches of 
heave. 
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DISCUSSION OF RESULTS 

The procedures to predict vertical differential soil movement were developed based 
on the breakthroughs mentioned earlier. These procedures consider moisture effect cases 
that are common with light commercial and residential structures, such as: (1) bare soil at 
the surface, (2) grass at the surface, (3) tree at the surface, and (4) flower bed at the surface. 
Additionally, these procedures include calculating the effects of differential soil movement 
caused by the introduction of design effects, such as, vertical barriers and horizontal 
barriers. The methodology to handle these moisture effect cases and design effects are 
quite delineated in detail in the process flow of detailed algorithms, Appendix A. A 
graphical representation of the program design layout is illustrated in Appendix B. 

The algorithms in Appendix A were programmed in a spreadsheet software 
package. The author used this package to analyze the moisture effect cases illustrated in 
Appendix C which consists of the soil profile described in Table 1. Table 4 is a summary 
of vertical differential soil movements for the site described by Appendix C. This 
summary illustrates the worst cases of edge lift and center lift distortion modes for all 
moisture effect cases. Post equilibrium and post construction approaches were used to 
solve for all cases so that the worst possible case of vertical differential soil movement is 
identified. Additional analyses were performed using different sizes of vertical and 
horizontal barriers to define the effects that these design features have on soil movement. 
An engineer may use the maximum cases of edge lift and center lift caused by the 
combination of moisture effect cases and the desired design effects to design an adequate 
foundation. 

Table 4 is set up in chronological order beginning with the post construction 
differential soil movements and concluding with the post equilibrium soil movements. 

The post construction numbers in Table 4 were generated for June as the driest 
month of the year and for a slab built in January (the wettest month of the year). The slab 
is placed on “wet” soil and this soil is allowed to approach an equilibrium moisture content 
based on the soil characteristics and the post construction theory described earlier. The soil 
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at the edge of the slab continues to fluctuate with the seasonal weather changes. Applying 
this post construction scenario to the “1-7 Tree” case produces the maximum differential 
soil movement compared to all other cases which cause a center lift distortion mode. The 
differential soil movement in this case is 4.32 inches and occurs the first August after 
construction. Note that a reduction to 2.44 inches of differential soil movement is made 
possible by introducing a four foot vertical barrier at the edge of the slab where the tree is 
located. 



Table 4. Summary of Vertical Differential Soil Movements for Appendix C 

Post Construction 

Driest Month of the June 

year 



Construction 



Worst Case Center Lift Distortion Mode 

Moisture Effect Case/Desiqn Scenario Differential Soil Movement Number of 

(units in inches) Months after 

Const. Worst 
Case Occurs 

1-4 Grass 2.34 6 

1-5 Grass with 4 ft Horizontal Barrier 1.19 6 

1-6 Grass with 4 ft Vertical Barrier 0.48 7 

1-7 Tree 4.32 7 

1 -8 T ree with 4 ft Vertical Barrier 2.44 7 

1-9 Flower Bed 3.30 6 

Worst Case Edge Lift Distortion Mode 

Moisture Effect Case/Desiqn Scenario Differential Soil Movement Number of 

(units in inches) Months after 

Const, Worst 
Case Occurs 

1-4 Grass 0.43 11 

1-5 Grass with 4 ft Horizontal Barrier 0.31 11 

1-6 Grass with 4 ft Vertical Barrier 0.23 14 

1-7 Tree 0.43 11 

1-8 Tree with 4 ft Vertical Barrier 0.23 14 

1-9 Flower Bed 0.92 12 

1-11 Flower Bed with 4 ft Vertical Barrier 0.20 12 
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Table 4. Continued 



Post Equilibrium 

Worst Case Differential Soil Movement 



Moisture Effect Case/Design Scenario 


Edqe Heave (in) 


Edqe Shrink 
(in) 


1-4 Grass 


1.65 


1.03 


1-5 Grass with 4 ft Horizontal Barrier 


0.47 


0.56 


1-6 Grass with 4 ft Vertical Barrier 


0.07 


0.09 


1-7 Tree 


1.64 


3.67 


1-8 Tree with 4 ft Vertical Barrier 


0.07 


2.36 


1-9 Flower Bed 


2.89 


1.86 


1-11 Flower Bed with 4 ft Vertical Barrier 


0.68 


0.55 


1-12 Measured Suction Profile (Point Water Source) 


4.40 


NA 


2-12 Measured Suction Profile (Point Water Source) 


4.40 


NA 


3-12 Measured Suction Profile (Tree Removal) 


1.62 


NA 



An inspection of Table 4 shows that the differential soil movement is greatest when 
there are no vertical or horizontal barriers in place. In the “1-4 Grass” case, the edge 
heave is 1.65 inches and edge shrink is 1.03 inches. When a four foot horizontal barrier is 
added the edge heave is reduced to 0.47 inches and the edge shrink is reduced to 0.56 
inches, as described in the Table 4 case “1-5 Grass with 4 ft Horizontal Barrier”. A further 
reduction in soil movement is obtained by adding a four foot vertical barrier as shown in 
Table 4 case “1-6 Grass with 4 ft Vertical Barrier”. This latter case reduces the edge heave 
to 0.07 inches and the edge shrink to 0.09 inches. 

As shown in Table 4, the worst case edge heave occurs in the post equilibrium case 
of the measured suction profiles produced at the locations of the leaking pipes. The 
leaking pipes produced a “wet” spot that caused the soil to heave 4.40 inches at those two 
locations. This type of heave at the edge caused the worst case edge lift distortion mode. 
Refer to Appendix B-130 for a conceptual illustration of the post construction case of 
vertical differential soil movement that can be predicted by applying the procedures 
contained herein. These procedures for analyzing vertical differential soil movements 
should be applied to a soil profile to enable an engineer to maximize resources ensuring 
that the best foundation is designed. 
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SUMMARY AND SUGGESTIONS FOR RESEARCH AND DEVELOPMENT 

The theory to predict differential soil movement of expansive soils has evolved 
since the 1970’s. The engineering profession is at a juncture today in which engineers 
have the technology to apply this information to maximize resources and produce high 
performing ground structures. The evolution of this theory began in the early 1970’s with 
the advent of the volume strain theory for expansive soils which was developed by Lytton 
(1973). In the decade of the 1970’s, great strides were made by Pearring (1968), Holt 
(1969), and Mojekwu (1979) to enable engineers to identify predominant clays. In 1981, 
McKeen used these newfound procedures to identify predominant clays and developed a 
correlation to link groups of clays to suction compression indices. The suction 
compression index is used in the soil volume strain equations developed by Lytton to 
predict vertical differential soil movement for particular types of expansive soils. In the 
1 980’s, although the volume strain theory was well based, its application was limited by 
the lack of knowledge and methodology to predict the equilibrium suction and soil suction 
profiles for a particular soil profile and location. This is when work by Mitchell and later 
Gay opened the door for major advancements in this area of expansive soils. In 1980, 
Mitchell developed and applied simple mathematical methods for predicting soil suction 
profiles. Then in 1994, Gay developed a way to predict the mean volumetric water content 
for soils dependent upon the location and climatic conditions. Jayatilaka, Gay, Lytton, and 
Wray (1992) then developed and calibrated to field data a finite element method for 
determining the edge moisture variation distance. They also devised methods to predict a 
number of unsaturated soil properties using simple index properties. 

Based on the wealth of information and breakthroughs in the areas of expansive 
soils listed above, Lytton and the author made five critical advancements necessary to 
effectively and easily apply this theory: 

1 . The means to calculate the equilibrium suction for a particular soil profile and location. 

2. The means to calculate depth to constant suction for a soil profile with multiple layers. 

3. The equation to calculate horizontal velocity of water flow in unsaturated soils. 
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4. A method to consider the transient case for calculating vertical differential soil 
movement for expansive soils immediately after construction and before the soil under 
the center of the slab has reached an equilibrium moisture content. 

5. Procedures to apply all of the theory discussed herein to predict differential soil 
movement for expansive soils considering all of the moisture effect cases and design 
effects that have been mentioned. 

The theory discussed herein represents physical behavioral characteristics of 
expansive and unsaturated soils and is well documented. Those contributing to this 
vertical differential soil movement theory have approached the problem with a 
theoretically sound yet practical mindset. Additionally, this expansive soil theory has been 
developed in a way to enable practicing engineers to run simple, common and economical 
geotechnical engineering lab tests to adequately define all necessary soils profile 
information. This soil profile information is applied in the vertical differential soil 
movement theory and should be used to make sound geotechnical engineering designs. 

Further development needs to be made in the areas of enhancements to the present 
programs, materials characterization testing, in situ testing, and foundation structural 
analysis models for stiffened plates and drilled shafts in expansive soils. The program 
enhancements should include the choice to use the Gardner, as well as the Mitchell, 
unsaturated permeability options. Materials characterization testing developments that are 
needed are more convenient methods of measuring the volume change, diffusion, and soil- 
water characteristics of expansive soils directly. In situ testing developments that are 
needed include methods of measuring suction profiles on site and preferably down hole; 
lateral earth pressure coefficient; and of detecting the presence of roots and voids, either 
air-filled or water-conducting; in the field and non-destructively. Stiffened plate models 
which incorporate non-rectangular geometry need to be coupled with this program or its 
enhancements to determine design quantities of moment, shear, and deflection. Drilled 
shaft models with interface elements need to be coupled with the enhanced program to 
determine the design quantities for both shrinking and expanding soils. 



78 



REFERENCES 

Castleberry, J. P. II, (1974). “An engineering geology analysis of home foundations on 
expansive clays.” M.S thesis, Dept.of Geology, Texas A&M University, College 
Station. 

Gardner, W. R., (1958). “Some steady state solutions of the unsaturated moisture flow 
equation with application to evaporation from a water table.” Soil Science, 85(4), 
223-232. 

Gay, D. A., (1994). “Development of a predictive model for pavement roughness on 
expansive clay.” PhD dissertation, Dept, of Civ. Engrg., Texas A&M University, 
College Station. 

Geostructural Tool Kit, Inc., (1996). PTISLAB WIN 1.0. Geostructural Tool Kit, Inc., 
Austin. 

Holland, J. E. and Lawrance, C. E. (1980). “Seasonal heave of Australian clay soils.” 
Proc., Fifth Int.Conf. on Expansive Soils. ASCE, Denver, 302-321. 

Holt, J. H., (1969). “A study of physico-chemical, mineralogical, and engineering index 
properties of fine-grained soils in relation to their expansive characteristics.” PhD 
dissertation, Dept, of Civ. Engrg., Texas A&M University, College Station, Texas. 

Jayatilaka, R., Gay, D. A., Lytton, R. L., and Wray, W. K. (1992). “Effectiveness of 
controlling pavement roughness due to expansive clays with vertical moisture 
barriers. ” Research Report No. 1165-2F, Texas Transportation Institute, College 
Station. 

Jones, D. E. and Holtz, W. G. (1973). “Expansive soils - the hidden disaster”, Civil 
Engineering , ASCE, August, 49-51. 

Juarez-Badillo, E. (1975). “Constitutive relationships for soils.” Proc., Symposium on 
Recent Developments in the Analysis of Soil Behavior and their Applications to 
Geotechnical Structures, The University of New South Wales, Kensington, N.S.W., 
Australia, 231-257. 

Lambe, T. W. and Whitman, R. V. (1969). Soil mechanics. John Wiley & Sons, New 



79 



York, New York. 

Lytton, R. L. and Woodburn, J. A., (1973). “Design and performance of mat 
foundations on expansive clay. ’’ Proc., Third Int. Conf. on Expansive Soils. 
American Society of Civil Engineers, Haifa, Israel, 301-307. 

Lytton, R. L., (1977). “The characterization of expansive soils in engineering.” 
Presentation at the Symposium on Water Movement and Equilibrium in Swelling 
Soils, American Geophysical Union, San Francisco. 

Lytton, R. L., Pufahl, D. E., Michalak, C. H., Liang, H. S., and Dempsey, B. J., (1989). 
“An integrated model of the climatic effects on pavements.” Rep. FHWA-RD-90- 
033 to Federal Highway Admin., Civ. Engrg., Texas A&M Univ., College Station. 

Lytton, R. L., (1994). “Prediction of movement in expansive clays. ” Geotechnical 
Special Publication No. 40, Yeung, A.T., and Felio, G.Y., eds., ASCE, New York, 
Vol. 2, 1827-1845. 

Lytton, R. L., (1995). “Foundations and pavements on unsaturated soils.” Proc., First 
Int. Conf. on Unsaturated Soils, E.E. Alonso and P. Delage, eds., Vol. 3, 1201-1220. 

Lytton, R. L., (1997). “Engineering structures in expansive soils.” Proc., Third 

Brazilian Symposium on Unsaturated Soils, De Campos, T.M.P., and Vargas, E.A. 
Jr., eds., Rio De Janeiro, Brazil, Vol. 2, 550-562. 

Mathewson, C. C., Castleberry, J. P. II, and Lytton, R.L., (1975). “Analysis and 
modeling of the performance of home foundations on expansive soils in central 
Texas.” Bulletin of the Association of Engineering Geologists, Vol. 12, No. 4, 275- 
302. 

Mathewson, C. C., Dobson, B. M., Dyke, L. D., and Lytton, R. L., (1980). “System 
interaction of expansive soils with light foundations.” Bulletin of the Association of 
Engineering Geologists, Vol. 17, No. 2, 55-94. 

McKeen, R. G. (1981). “Design of airport pavements on expansive soils.” U.S. 

Department of Transportation, Federal Aviation Administration, Washington D.C. 

McQueen, I. S. and Miller, R. F., (1968). “Calibration and evaluation of a wide-range 
gravimetric method for measuring moisture stress.” Soil Science, Vol. 106, No. 3, 



80 



225-231. 

Mitchell, P. W. (1980). The structural analysis of footings on expansive soil 2 n ^ 
edition , K.W.G. Smith and Associates, Newton, South Australia. 

Mojekwu, E. C. (1979). “A simplified method for identifying the predominant clay 
mineral in soil.” M.S thesis, Dept, of Civ. Engrg., Texas Tech University, Lubbock. 

Nieber, J. L., (1981). “Simulation of infiltration into cracked soils.” Presented at the 
1981 Summer Meeting, American Society of Agricultural Engineers, Orlando. 

Pearring, J. R, (1968). “A study of basic mineralogical, physical-chemical, and 

engineering index properties of laterite soils.” PhD dissertation, Dept, of Civ. Engrg., 
Texas A&M University, College Station. 

Post-Tensioning Institute, (1980). Design and construction of post-tensioned slabs-on- 
ground, 1 st edition. Post Tensioning Institute, Phoenix. 

Post-Tensioning Institute, (1996). Design and construction of post-tensioned slabs-on- 
ground, 2 n> ^ edition. Post Tensioning Institute, Phoenix. 

Skempton, A. W., (1955). “The colloidal activity of clays. ” Proc., Third Int. Conf. On 
Soil Mechanics and Foundation Engineering, Vol. 1, Switzerland, 57. 

Spotts, J. W., (1974). “The role of water in gilgai formation.” PhD dissertation, Texas 
A&M University, College Station. 

Thomthwaite, C. W., (1948). “An approach toward a rational classification of 
climate.” The Geographical Review 38, 55-94. 

Viayvergiya, V. N., and Ghazzaly, O. I., (1973). “Prediction of swelling potential for 
natural clays.” Proc., Third Int. Conf. on Expansive Soils, Haifa, Israel, Vol. 1, 227- 
236. 

Wiggins, J. H., (1976). “Natural hazards, an unexpected building loss assessment.” 
Technical Report No. 1246, J. H. Wiggins Co., Redondo Beach, California. 

Wray, W. K. (1978). “Development of a design procedure for residential and light 
commercial slabs-on-ground constructed on expansive clays.” PhD dissertation, 
Texas A&M University, College Station. 



81 



APPENDIX A 

PROCEDURES TO CALCULATE 
VERTICAL DIFFERENTIAL SOIL MOVEMENT 
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Start of Program 
and 





Note to Programmer: Most of the 
input should draw from 
databases or spreadsheets in a 
"Combo Box" format to provide a 
user friendly layout. 










\ * 


Project Information Input 

1. Input Project Name 

2. Input Project Date 

3. Input Project Engineer 

4. Input Site Location Address 


... 


/ * 






/Note to Programmer: This 
information should be stored in a 
database file of "Projects" in a 
spreadsheet type format. This 
"Projects" file could then be used 
> to organize all projects analyzed 
\by this software package. 
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Foundation Layout Input General Note 

Note to user:A three dimensional rectangular (x,y,z) 
coordinate system is used to define foundation layout. The Y 
direction is defined by an angle wrt true north. The z direction is 
depth of soil (increasing positive below surface). Every part of the 
\ foundation structure must be in the positive (x,y) plane and the slab 
shape must be rectangular. The origin shall be defined by latitude, 
longitude, and elevation. Refer to Appendix B-118, 125, and 126 
\ for graphical representation. 




Input the slab origin 
(typically the west- 
most and the south- 
most corner of the 
slab. 



Input the elevation of 
the origin. 



Input the y axis 
direction 



Input the Local TMI 
TMI = 





j 


Input the number of 
Local weather cycles/yi 

n = 


y 


f - — 



Input the information 
for the corners of the 
slab. 




Specific Notes 



Inputs for latitude and longitude 
are intended for data collection 
purposes and to link the site to site 
specific data, such as; Thornthwaite 
Moisture Input data, and edge 
moisture variation data. 

Refer teethe example below. 



The origin elevation is simply used as 
a reference for 3D computations. 

Refer to the example below. 



This hasno bearing on calculations. 
This simply allows the user to establish 
building location/alignment. The value 
should be entered clockwise from true 
north. Refer to the example below. 



input Origin; Latitude | 29°15’ N 



Longitude |96^2?W 

'■'■‘•5 •• 'V 

Origin Elevation I 100 00 (ft-) 

■ I I 

Angle Clokwise from 
True North 



Input Direction of “Y* Axis | 4°30‘ 



Example: Cut-out 1 from Appendix B-118 



Specific Note 



All cells in the "corner information" 
spreadsheet are critical. The user 
must not be allowed to continue without 
first entering this information. 
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Soil Layer Property 
Input 




No 




Calculate additional 
soil property quantities 
required to make 
additional 
computations. 




◄ — ( from ) 
\pg88/ 



Give user the option to view 
spreadsheet of other calculated 
material properties used in volume 
strain, suction and water flow 
velocity calculations. Refer to 
Appendix B-126. 
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Moisture Effect Input 




Example: Cut-out 2 from Appendix B-121 



Program Design Layout 



INPUT j OUTPUT | 



ThmtitimiU MU# | 0*m | MUnx 



J 



20 4p J~ 



leaking Pip. 
pF 2.0 



Vertical 
p arris r 

length pf ih* East edge el 



60 60 100 



Example: Cut-out 3 from Appendix B-125 



Start Location of Effect End Location of Effect 



X Coord 



1 



Y Coord. X Coord. Y Coord. 



T 



Example: Cut-out 4 from Appendix B-121 



As the Moisture Effect Case types 
are selected a database 
spreadsheet row is generated anc 
a counter is initiated for each case 



Example 



For example, the first "flowerbed" moisture 
effect case selected would be given the 
"case number" name (1-6) (and the second 
would be (2-6) and so on]. A row would be 
established to contain all pertinent cells of 
information required for computations made 
by this software package for that particular 
case. Refer to cut-out 7 on the next page. 



Prompt user to input a "start" and "endp 
location for all moisture effect cases. 



Example 




For example, the "start" and "end" points for 
the flower bed are (40,0) and (0,50) 
respectively. If the moisture effect is a point 
source then the two inputs will be the same. 
Refer to the cut-out 3 and cut-out 4 for 
location of flowerbed on grid and 
example of coordinate inputs respectively. 
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Moisture Effect Input 
Continued 



Select & Input Molettxre Effect 
Case 



Measured Sue tin rt Profile 
Input 



Select Type of Moisture Effect Case 
Start Location of Effect End Location of Effect 




| X Coord, | 


| Y Coord. | X Coord. j| Y Coord. 


1 40 I 


1 o 


1 0 II 50 


jlinttort k*»t«d «! edf* ed founckt* 


a£] Irput Edge Moisture Variation 


Depth of Eflfcet 
(ft.) 




| Edge Lift || Center Lift j 


1 4.0 j| 


II _ 

* 


1 1 .Jill IS 



IT Ho", ditunci from «4pt 

(W»p«U»» R**b«r (kfoaitju noin vn «U*ct «df* to U wder tt» fcxaditfei ) 



C 7 A 

( from ! 
Vpg 85 



1 Bln soil »t iurt»c« 

2 Bln roilvuitbboriioMilbiniir 

3 Eire »il u4th verticil tenter 
4. On u at turtle* 

5 On u vault luruoiti] barrier 

6 Gn u unth vtrucil binder 

I Tn« near foundation *dr* 

9 Tn« uutli verticil burner 

9 . Flowerbed near foundation «dje 

10. Flower b«d with koritontilbiniar 

I I Flower b«d uitX ee rtieil bi m*r 
12 Measured sun son profile i 



Perform the following 
function for moisture effect 
cases 7 thru 1 1 . 




Example: Cut-out 5 from Appendix B-121 



For example, the depth of the flower bed in 
this case is 4 feet. Refer to cut-out 5 shown 
above. 



Refer to cut-out 6 for example of edge 
location input for the flower bed. 



Example 



Example 



Prompt user to input the depth of the 
moisture effect zone. 
d„ 



Is effect located it ed ge of found* lion? T 




distance from edge. 

(Negative number defines the 
moisture effect edge to be under the 
foundation.) 






> 




Prompt user to enter 
location wrt slab. 

(-) number indicates 
the zone is under the 
slab, (+) number 
indicates the zone is 
away from the slab. 



When yes is selected 
this indicates the effeci 
is at the edge of the 
slab. 



Example: Cut-out 6 from Appendix B-121 



Provide the user opportunity to 
review spreadsheet of moisture 
effect input data. Refer to cut- 
out 7. 




’K 

J t lx Coord 


“1 
Y Coord 


End | 

I Coord | YCoont j 


Barrier 

Eire 

<«-) 


fejrtb Edge I 
MftHUgaUotU' 


1 is 1 « 


56 


65 f 120 | 


1 .. 4 .. 


1 «J ...» I| 


P7 I 


0 | 


75 


o | no j 




1 20 J .7.5 J| 


« 1 


40 


o 


.03 50 ;j 




1 4 1 “ l| 


1 H. L 65 


120 


Vi <n 


4 


1 1 .0 1 


P* | 


65 


100 


— r — 1 — n — i 


1 2-12 | 


65 


20 


r r 1 1 1 1 


| 3 45 


60 


SI i 


1 i 


r~i — 1 




6 A 

go to * 
pg 85/ 



Example: Cut-out 7 from Appendix B-121 
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Measured Suction 
Profile Input 




88 
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Calculations for 
Post Equilibrium 
Moisture Effect 
Cases 



Calculations for 
Post Construction 
Moisture Effect 
Cases 




— Select 2 


f 


Input Construction 
Time Constraintsand 
establish Time 
Baselines 












r 


Calculate the Suction 
Profiles for all Moisture 
Effect Cases. 









Calculate the 
Differential Soil 
Movement for all 
Moisture Effect Cases. 
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Equilibrium Suction 
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Depth to Constant 
Suction 




The algorithms used to compute the 
depth to constant suction may vary 
depending on the moisture effect 
case being considered. 
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Depth to Constant 
Suction 
For Bare Soil 



Depth to Constant 
Suction 
For Grass @ 
Surface 
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Depth to Constant 
Suction 

For a Measured 
Profile 





Define Variables 

s = layer number considered given an initial value of zero. s=0 
L depth to layer s (retrieved from pg 84) 
z s = depth of increment s (in cm), 

z m = depth to constant suction given an initial value of zero. ^=0, 
h c = the c lh suction value for a sample (retrieved from pg 87) 

N= Number of suction measurements per sample (retrieved from pg 8 
D n = depth to the N th suction measurement (retrieved from pg 87), 
h rn (D N )=equilibrium suction at depth IJJ (retrieved from pg 96) 

U 0 = absolute value of (h N - h m (D N )) 

Z o = °N 
^Mitl = 00 













N 

3 

II 




z m =z m value for the surface 
conditions most similar to 
the surface conditions for 
the measured profile but 
not less than ^ for bare 
soil. 




i 

Yes 1 




i 


i 




Depth to Moisture Active 
Zone 

Z m “ Z m + ^final 
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Calculate Suction 
Profiles 
Menu Page 





The algorithm used to calculate a 
suction profile wrt depth varies 
depending on the moisture effect 
case and surface conditions. 



Go to the Propersuction profile 
algorithm depending on which moisture 
effect case is being considered. 
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97 



Dry Limit Suction 






98 



Dry Limit Suction 
Profile 

For Tree and Flower 




99 



Wet Limit Suction 




Wet Limit Suction 
Profile 

For Flower Bed (a) 





Yes 
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Suction Profile 
For a Measured 
Profile 




102 



Calculate Soil 
Volume change 
Menu Page 




The algorithm used to calculate soil 
volume change varies depending 
on the moisture effect case and 
surface conditions. 
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Soil Volume Change 
for Bare Soil, Grass, Tree 
and Flower bed at the 





104 



Soil Volume Change 
for Bare Soil, Grass, Tree and 
Flower Bed at Surface with 
Horizontal Barrier 




105 




36 (a) N 
go to 
pg 106/ 
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Soil Volume Change 
for Measured Suction Profiles 
(Continued) 
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Post Construction Time 
Constraint Input and 
Calculations 



Definition of Variables 
t dry = Driest month of the year. 
t const = The month the slab was constructed, 
t^ The month the soil volume change is sought. 

^ const/dry = Time elapsed between the driest month and 
and the month of construction, 
tyc/const = Time e * a P secl between the months of 

construction and the month the soil volume 
change is sought. 

tyc/dry = Time e,a P sed between the driest month and th 3 
month the soil volume change is sought. 

test1= Dummy variable. 

test2= Dummy variable. 




Input the driest month of the year t dry = _ 

Input the month of construction t const = 

Input the month soil volume change is sought t vc = 



Set text values equal to 
numberical values fort dry) 

^const* and *vc 



"Jan"=0 


"Jul"=6 


"Feb"=1 


"Aug'-T 


"Mar"=2 


"Sep'-B 


"Apr"=3 


"Ocr=9 


"May"=4 


"Nov"=10 


"Jun"=5 


"Dec'-I 1 




test2 = 12 



test2 = 0 






testl = Ust - t dry 





r 


■ ‘estl = t„ c - t corKt 
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Calculate Post 
Construction 
Suction Profiles 
Menu Page 
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Retrieve and 
Intaialize Variables 
for Suction Profile 
Algorithms 




Bare Soil at Surface. 
Retrieve Values from previous 
process flow sheets: 

tvc/dry ( SeC > fr0m P9 107 ' 

Z m fr ° m P9 92, 
h m from pg 90, 
hJZj) from pg 96, 
n from pg 83, 
a(Zj) from pg 88, 
h dry from pg 97, 

h m a x< z i) from P9 97 
h w , from pg 99, 

M b from pg 84, 

L s from pg 84, 

^ _ ^vc/dry* 





42 



Grass at Surface. 
Retrieve Values from previous 
process flow sheets: 

tvc/dry ( SeC > fr0m Pg 107 > 

z m from pg 92, 
h m from pg 90, 
h m (Zj) from pg 96, 
n from pg 83, 
a(Zj) from pg 88, 
h dry from pg 97, 

from pg 97 * 

h wl from pg 99, 

M b from pg 84, 

L from pg 84, 



Tree at Surface. 
Retrieve Values from previous 
process flow sheets: 

fyc/dry ( SeC ) fr0m Pg 107 > 

z m from pg 93, 
h m from pg 90, 
h m (z.) from pg 96, 
n from pg 83, 
a(z.) from pg 88, 
h from pg 98, 

h max( z i) from P9 98 . 
h wl from pg 99, 

KninW from "• 

M b from pg 84, 

L s from pg 84, 





Under Slab Profile at time of 
construction. 

Retrieve Values from previous 
process flow sheets: 

tconst/dry ( SeC > fr0m P9 107 ' 

t = t 

const/dry 

Use the time variable listed 
above (t=t const/dry ) in the place of 
time variable (t^^) and 
perform suction profile 
calculations depending on the 
moisture effect conditions that 
existed at the surface just prior t{> 
placing the slab. 



Flower Bed at Surface. 

Retrieve Values from previous 
process flow sheets: 

t vc ( dn,< SeC ) fr0m p 9 107 - 
z m from pg 93, 
h m from pg 90, 
h m (Zj) from pg 96, 
n from pg 83, 
a(z.) from pg 88, 
h dry fr0m pg 98 > 
bmax^i) fr0m Pg 98 ’ 
h wl from pg 100, 

M b from pg 84, 

L s from pg 84, 

^ — ^vc/dry* 
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Suction Profile 
Algorithms 




Ill 





1 12 



Post Construction 
Soil Volume Change 
Menu Page 






The algorithm used to calculate soil 
volume change varies depending 
on the moisture effect case and 
surface conditions. 



Go to the Propersoil volume change 
algorithm depending on which moisture 
effect case is considered. 



Calculate 
Volume Change 
for the following 



cases (if the case exists): \ 

1 . Bare soil at surface, 

2. Bare soil with horizontal barrier, \ 


1,2,3 


3. Bare soil with vertical barrier, \ 

4. Grass at surface, \ 


5. Grass with horizontal barrier, 

6. Grass with vertical barrier, 


\ 4,5,6' 


7. Tree at surface, 

8. Tree with vertical barrier, / 


/ X. 7,8,12 


9. Flower bed at surface, / 


9,10 


10. Flower bed with horizontal barrier, / 

1 1 . Flower bed with vertical barrier, / 

12. Moisture effect some distance from/ 

the barrier/slab edge, / 

13. All quantities calculated. / 


11,12 







13 
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Retrive Values for a 
Particular Case of Post 
Construction Soil Volume 
Change Calculations 
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Post construction Soil 
Volume Change 
Calculations With or 
Without Barriers 
(Part I of II) 

00* No - is ^ <=z " ? ' 



/ 53 \ 
from 
P9113/ 



i = i + 1 



[ 

Yes 



Convert all Suction in pF 
to centimeters. 
h(cm) = 1 o h(pF> 



▼ 



Az^ = Az cd(;e + increment 






| <: ^ S J 1 ^'^dr^ > ^slab^'^dry^^ 






Yes 

I 

* 

K = 33 
f = .5 

e mc =Center fr0m P9 88 



No 



k 0 = .67 
f = .8 

e ce = edge lift from pg 88 



if] 


1 =0 
i 




(AV'j ( AV^ ( AV'l 

IitJ.-I-vJ/IitJ. 



15 Z i >=d vbar ? 




Calculate Horizontal Velocity 

.. , . _ Ko( Zi ).|h 0 | fh slab ( Zi , W 

( e mc/c + d sourcc) V h ( Z i ’ 1 vc/dry ) . 



Calculate Suction at Foundation edge 



^edge( Z i » ^ vc/dry) ^slab( Z i ’ ^ vc/dry ) * ^ 



Mzi)k| 




^F^ 2 i) “ ^edge( Z i’ ^vc/dry) 



h F ( Z j) “ ^stab^ Z i’^/c/dr^ 



(^) =-Yh(Zi)>ogi.i 



^l( Z |) ~ ^slab^ Z i’^vc/dry^ 






▼ 




M z i) = 


1 + 2-ko" 

3 


’ z i 



CT F (Z,) = 0 
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Post construction Soil 
Volume Change 
Calculations With or 
Without Barriers 
(Part II of II) 



Az siab = Az siab + increment- 



I 



(IHf). 










j = j 


+ 1 




Convert all Suction in pF 
to centimeters. 
h(cm) = 10 h<|,F) 


Zj = z j _ 1 + increment 





^const/dry) > ^slab^j’^vc/dn 


yr^ 


No 




Yes 








CO 

CO 

II 

c 




k 0 = .67 


f = .5 




f = .8 



Yes 




I'll 2 ,) ^slab( Z i’tconst/dcy) 



h F< Z i> = h slab( z i' W 



-r h ( Z ,) ,0 g„(^j 









( AV'j 

^ vJ„ 


-($ 


<T 


r^>YeS| 


r 




■1% 


) 

a 
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APPENDIX B 

SOFTWARE PROGRAM DESIGN LAYOUT 
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Program Design Layout 



VOLFLO WIN 1.0 



BH1Q 




INPUT 



GRAPHIC 

klSEIAlL- 



OUTPUT 






1 18 



Program Design Layout 





1 19 



Program Design Layout 
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Program Design Layout 
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Program Design Layout 
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Program Design Layout 




INPUT 






GRAPHIC 

BISEIi&L 



INPUT 



OUTPUT 



Soil Layer Prop. Input 



Project Into. 



Moisture Eflfect Input 



Foundation 



Select & Input Moisture Effect 
Case 



Measured Suction Profile 
Input 



Input Measured 
Suction Profiles 



Input Number of 
Measurements per Profile 

Note: Equilibrium suction will be assumed as default values for the 
initial/ final if nothing is entered. 



Case 

Number 



Number 
of Suction 
Me as urements 
Per Sample 



Initial Final Surface 

Suction Suction Conditions 

Profile(pF) Profile(pF) (Vegetation/Bare) 



1-12 


5 


It.II Input 


bi 


Input | 


Vegetation 


242 


5 


^|| Default 


IH 


Default 1 


Vegetation 


342 

L 1 


5 


v-J Input 


a 


Input 


Vegetation 



M ?£ '# •$> A 4V & i% 
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INPUT 



Program Design Layout 



INPUT 


GRAPHIC 1 OUTPUT 1 




DISPLAY 1 | 



Project Into. | Foundation Layout | Moisture Effect Input j Soil Layer Prop. Input | 
ffect | 



Select & Input Moisture Effect 
Case 



Input Number of 
Measurements per ProBie 



Measured Suction Profile 
Input 

Input Measured 
Suction Profiles 



Case 

Number 


Suction 

Measurement 

Number 


Depth to 
Measurement 
(ft) 


Suction 

Measurement 

(pF) 


Initial/ 

Final 


1-12-F 


1 


0.0 


3.0 


Final 


1-12-F 


2 


3.0 


2.0 


Final 


1-1 2-F 


3 


12.0 


3.0 


Final 


1-12-F 


4 


16.0 


3.1 


Final 


I 1-12-F 


5 


20.0 


3.3 


Final 


2-1 2-F 


1 


0.0 


3.1 


Final 


2-1 2-F 


2 


3.0 


1.8 


Final 


I 2-1 2-F 


3 


12.0 


3.0 


Final 


| 2-1 2-F 


4 


16.0 


3.1 


Final 


2-1 2-F 


5 


20.0 


3.2 


Final 


3-12-1 


1 


0.0 


4.5 


Initial 


3-12-1 


2 


3.0 


4.5 


Initial 


3-12-1 


3 


12.0 


4.5 


Initi al 


I 3-12-1 


4 


16.0 


4.0 


Initial 


3-12-1 


5 


20.0 


3.3 


Initial 
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APPENDIX C 

SITE PLAN OF MOISTURE EFFECT CASES 
USED IN SAMPLE CALCULATIONS 
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APPENDIX D 

REQUEST FOR PERMISSION TO USE COPYRIGHTED MATERIAL 
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Donald D. Naiser, Jr. 

400 Marion Pugh #1905 
College Station, Texas 77840 



August 27, 1 997 



Permissions Department 
John Wiley & Sons, Inc. 

605 Third Avenue 

New York, NY 10158-0012 

REQUEST TO USE FIGURES COPIED FROM SOIL MECHANICS 1969 

I am a master of science candidate at Texas A&M University and I will be finalizing my 
thesis in the next two weeks. Our library sends our theses to University Microfilms Inc., for 
preparation of a microfilm copy of the document. As I am sure you are aware, UMI retains the 
right to sell both hard and microfilm copies of the document. 

My reason for writing is to request permission to use the figures listed below exactly as 
they appear in SOIL MECHANICS , by Lambe and Whitman, copyrighted by John Wiley & 
Sons, Inc., 1969. With your permission, I will use the following figures in my thesis: 

1 . Figure 4.7 “The structure of serpentine.” Found on page 47. 

2. Figure 4.8 “The structure of kaolinite.” Found on page 47. 

3. Figure 4.9 “The structure of pyrophyllite.” Found on page 47. 

4. Figure 4.10 “The structure of muscovite.” Found on page 47. 

5. Figure 4.1 1 “Sheet silicate minerals.” Found on page 48-49. 

6. Figure 5.9 “Soil particles with water and ions.” Found on page 55. 

One chapter in my thesis describes the magnitude of the problem of designing 
foundations built on expansive soils. The figures described above will be used to illustrate the 
molecular structure of fine grained soil particles. The thesis has been accepted and the author 
proofs indicate a 1997 publication date. The thesis title is PROCEDURES TO PREDICT 
VERTICAL DIFFERENTIAL SOIL MOVEMENT FOR EXPANSIVE SOILS. 

Please inform me of any necessary action to prevent future publication complications. I 
can be contacted by phone at home (409) 543-8412 or at work (409) 845-9919. Thank you very 
much for your effort in resolving these issues. 



Sincerely, 




Donald D. Naiser, Jr. 
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VITA 

Donald David Naiser, Jr. is a citizen and resident of the United States of 
America, whose permanent mailing address is, 

Rt. 3 Box 203-C 
El Campo, Texas 77437. 

He obtained his Bachelor of Science degree in Civil Engineering, from Texas 
A&M University in 1991. 

Since graduating from Texas A&M University, he has been commissioned as a 
United States Naval Officer as a member of the Civil Engineer Corps. His time served in 
the Navy has been exclusively with the SEABEES performing contingency construction 
training as well as constructing naval facilities. His tour highlights have been as the 
Officer In Charge of the Elevated Causeway System (ELCAS) in Amphibious 
Construction Battalion ONE and as the Officer In Charge of Construction Battalion Unit 
418. 

In 1996 he was accepted at Texas A&M University to pursue studies toward the 
degree of M.S. within the Geotechnical branch of Civil Engineering. 

In 1997 he became registered as a professional engineer in the state of 
Washington. 

“Blessed are the people whose leaders can look destiny in the eye without 
flinching but also without attempting to play God.”, Henry Kissinger. His career 
interests include achieving excellence through leadership in any Christian form. 
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